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ABSTRACT 


In  the  differential  equation 


h dy 

E A 

dx 


A(x,  e)y 


let  A(x,  e)  be  an  n x n matrix  valued  function  of  the  complex  variable  x 
and  the  parameter  e,  holomorphic  in  both  variables  for  |x|  < xQ,  0 < e < eQ. 
Let  y be  a vector  and  h a positive  integer. 


Fundamental  matrix  solutions  of  the  differential  equation  are 
constructed  which  involve  functions  that  can  be  represented  asymptotically 
by  series  in  powers  of  e.  The  domains  of  validity  of  these  expansions  grow, 
as  e tends  to  zero,  in  such  a way  that  their  distance  from  x = 0 tends  to 
zero  with  e.  Combined  with  stretching  transformations  of  the  form  x = £ep, 
this  result  is  a basic  tool  for  the  analysis  of  turning  points  of  unrestricted 
complexity. 

Equivalent  results  are  contained  in  papers  by  M.  Iwano  (Items  [6], 

[ 7 ] , [ 8 ] of  the  bibliography  at  the  end  of  this  paper. ) The  method  of  the 
present  article  has  many  points  of  similarity  with  Iwano 's  work,  but  it  is 
substantially  simpler  and  shorter.  This  is  achieved  by  a systematic  use  of 
a somewhat  generalized  concept  of  asymptotic  series  in  powers  of  e for 
functions  of  x and  e,  by  drawing  on  results  by  H.  Turrittin  and 
V.  I.  Arnold,  and  by  avoiding  the  preparatory  reduction  of  the  problem  to  a 
block-triangular  form,  necessary  in  Iwano's  presentation. 
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OUTER  SOLUTIONS  FOR  GENERAL  LINEAR  TURNING  POINT  PROBLEMS 


Wolfgang  Wasow 


§ 1 . Introduction. 

There  is  a vast  body  of  literature  on  the  properties  of  linear  analytic 
differential  equations  near  a pole  with  respect  to  a parameter.  For 
general  investigations  such  differential  equations  are  usually  written 
in  the  form 

(1.1)  eh  dJ  " A(x,  e)y, 

where  A(x,  e)  is  an  n xn  matrix- valued  function,  y an  n-dimensional 
column  vector,  e a small  parameter  and  h a positive  integer. 

The  main  difficulty  in  these  theories  stems  from  the  fact  that  the 
algebraically  relevant  properties  of  the  coefficient  matrix  A(x,  e)  may 
change  discontinuously,  even  at  points  where  it  is  holomorphic.  Some 
restrictions  on  A(x,  e)  have  therefore  been  imposed  in  most  papers  in 
this  field,  whether  they  deal  with  local  or  global  questions. 

M.  Iwano,  in  three  connected  papers  [fa],  [7],  [8],  has  developed 
a theory  that  is  truly  general,  in  that  for  any  problem  of  the  form  f 1 . 1)  a 
method  is  described  by  which,  for  any  point  x = xQ,  a set  of  asymptotic 
series  for  certain  fundamental  solutions  car.  be  calculated  so  that  the 
domains  of  validity  of  these  expansions  together  cover  a whole  neighborhood 
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of  xQ.  The  great  length  and  complexity  of  these  investigations  makes 
it  desirable  to  simplify  and  shorten  the  presentation,  so  as  to  make  it 
more  easily  accessible.  This  is  the  aim  of  the  present  paper.  Its 
arguments  have  many  points  of  similarity  witn  Iwano's,  but  there  are 
also  some  important  differences.  In  particular,  I permit  from  the 
beginning  matrices  that  are  not  block-triangular,  a restriction  which  Iwano 
removes  only  in  the  Appendix  of  his  third  paper.  The  widened  concept  of 
asymptotic  power  series  introduced  in  section  2 is  a helpful  tool  which 
shortens  the  subsequent  analysis.  In  addition,  the  presentation  is 
facilitated  by  employing  certain  results  of  Turrittin  [ 17 1 and  of  Arnold  [ 1] . 

Near  points  of  the  complex  x-plane  where  ail  eigenvalues  of  A(x,  0) 
are  distinct,  a complete  solution  of  (1.1)  in  terms  of  asymptotic  series^ 
in  powers  of  e has  been  known  since  early  in  this  century.  Its  descrip- 
tion can  be  found  in  several  textbooks,  e . g . , [ 2 ) or  ( 18 ) . 

The  existing  knowledge  near  points  where  A(x,0)  has  multiple 
eigenvalues  is  still  incomplete.  If  A(x,0)  has  no  eigenvalues  which 
are  identically  equal,  there  exists  a set  of  isolated  points  in  the  x-plane 
at  which  multiple  eigenvalues  occur.  They  are  often  called  "turning  points". 
The  asymptotic  nature  of  the  solutions  near  such  points  can  be  extra- 
ordinarily complicated  (see,  e.g.,  [9],  [ 10  J , [ 11 ))  and  involves  as 
yet  unsolved  connection  problems.  If  A(x,  0)  possesses  some  eigen- 
values that  are  identically  equal,  it  is  not  even  clear  what  terminology 
should  be  adopted.  Is  it  then  reasonable  to  call  every  point  in  the  domain 
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of  A(x,  0)  a turning  point?  One  might  be  tempted  to  apply  that  term 
only  to  points,  where  some  eigenvalues  are  equal  without  being  identically 
equal.  However,  as  Hukuhara  [ S ] and  Turrittin  [17]  have  shown,  the 
asymptotic  analysis  near  points  where  all  multiple  roots  are  identically 
equal  often  leads  unavoidably  to  problems  that  have  there  a turning 
point  in  the  former  sense. 

I shall,  therefore,  not  attempt  to  offer  a mathematically  precise 
definition  of  the  concept  of  "turning  point",  but  I will  use  these  words 
in  a descriptive  sense  as  having  to  do  with  difficulties  in  the  asymptotic 
theory  arising  near  certain  points  from  abrupt  changes  in  the  algebraic 
structure  of  relevant  matrices  at  that  point. 

The  most  generally  applicable  method  of  approach  to  turning  point 
problems  is  based  on  "stretching"  and  "matching".  If  the  point  in 
question  is  at  x - 0,  a transformation  of  the  form  x = £eK,  k > 0, 
is  fittingly  called  a stretching.  For  properly  chosen  values  of  k the 
"stretched"  differential  equation  can  often  be  solved  or,  at  least, 
brought  closer  to  a solution,  in  bounded  domains  of  the  £,-plane.  As  the 
image  in  the  x-plane  of  such  a domain  shrinks  to  zero  with  e,  a 
successful  "matching"  requires  that,  in  addition,  asymptotic  solutions 
are  available  in  regions  which,  while  possibly  bounded  away  from  the 
turning  point,  expand  towards  it,  as  e tends  to  zero.  The  two  types 
of  solutions  are  often  called  "inner"  and  "outer"  solutions. 


- 


As  Iwano  and  Sibuya  have  shown  ( | (■  ] , [ 7 ] , | 8 ) , [ 10  ] ) the  stretched 


problem  may  again  require  a composite  analysis,  involving  series  of 
outer1'  and  "inner"  type  with  respect  to  the  new,  stretched,  variable. 

The  ingenious  and  methodical  work  in  the  papers  cited  above  does  not 
yet  solve  the  problem  of  the  connection  of  the  several  fundamental  solu- 
tions with  known  asymptotic  series  valid  in  different  domains.  For 
partial  results  concerning  the  calculation  of  the  "connection  matrices' 
see  | 14  | and  [21].  As  far  as  the  calculation  of  each  of  these  solutions 
is  concerned,  all  the  features  appear  already  in  the  outer  solution  for  the 
original  differential  equation.  A few  further  remarks  on  this  point  can 
be  found  at  the  end  of  section  11. 

Because  of  the  length  of  the  arguments  it  is  desirable  to  begin  with 
a description  of  the  main  result  in  a language  that  is  not  excessively 
technical: 

The  asymptotic  solution  of  the  differential  equation  (1.1)  constructed 
in  this  paper  involves  asymptotic  series  in  powers  of  £ whose  coefficients 
are  not  necessarily  holomorphic  in  x at  the  point  x = 0 near  which 
the  outer  solutions  are  to  be  found.  These  coefficients  may  grow  as  fast 
as  some  negative  powers  of  x,  as  the  point  x = 0 is  approached. 
However,  the  exponents  of  these  negative  powers  do  not  increase  faster 
than  at  a linear  rate  with  the  order  of  the  term  in  the  series.  Thus,  if 
f (x)er  is  the  rth  term  of  sucii  a series,  there  exist  two  numbers,  p 
and  k,  independent  of  r such  that  fr(x)xP+r^K  remains  bounded  as 


-4- 


x - 0,  at  least  in  some  sector  of  the  x-plane.  Functions  of  x and  e 
that  possess  such  expansions  are  said  to  be  in  class  (7  . For  precise 

definitions  read  Definitions  2.1  and  9.1  below. 

Theorem  1.1  (Main  Theorem).  Let 

(1- 1)  eh  dx  = e^y:  ^ 51  °>  integer  , 

be  a system  of  n scalar  differential  equations  whose  coefficient  matrix 
is  holomorphic  in  lx  I < x , 0 < c £ eQ  and  possesses  an  asymptotic 
expansion 


00 

A(x,  e)  ~ Yi  e - 0+  , 

r=  0 

(1.2)  uniformly  in  i arg  x I < 6Q.  Then  (1.1)  has  a fundamental  matrix 
solution  Y(x,  e)  of  the  form 


(1.  3) 


Y(x,e)  = Y(x,e)eE  Q(x’ e) 


with  the  following  properties 

~ 1/m  1/1  ' t .....  „ . 

(i)  x = x , e = e ; m,  t are  positive  integers.  H is  a non- 


negative integer. 


(ii)  Y(x,  e)  and  Q(x,  e)  are  matrices  in  class  d in 
I arg  x I < 6^  < 6Q. 

(iii)  Q(x,  e)  is  a polynomial  in  e. 

(iv)  Q(x,  e)  is  diagonal. 

(v)  The  determinant  of  the  leading  term  YQ(x)  in  the  expansion  of 


Y(x,  e)  is  not  identically  zero. 


The  method  of  proof  is  completely  constructional.  This  does  not 
exclude  that  for  many  concrete  examples  shorter  procedures  may  be 
available.  Most  of  the  required  operations  are  rational.  There  occur, 
in  addition,  calculations  of  eigenvalues  of  holomorphic  matrix  functions 
and  quadratures  of  analytic  functions. 

The  by  now  fairly  common  notation  A = : B will  occasionally 
be  used  to  indicate  that  the  formula  represents  a definition  of  B.  The 
meaning  of  A :=  B is  analogous. 
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(uj)  There  is  a constant  k,  0 < k £ » such  that 


(2.  4) 


(2.5) 


11  r < r/x,  r = 0,1,...  . 


(iv)  Define  the  domain  f bv 

e 


£ = {x 

e 


-1  K 


0 e < lx  I < xQ,  I arg  x - 0Q  I < 6Q}  , 


where  tQ  >0,  xy  > 0,  6^  > 0,  and  0^  are  certain  constants.  Then 
there  is  for  every  N > 0 a constant  c^,  depending  on  N,  but  not 


on  x or  e,  such  that 

(2.  6) 


N 


i A(x,  e)  - £ A(x)er|<c  (|x|  1//Ke)N+1  , 
r=0  r 


for  x t f.  , 0 < e < c„. 

— e - 0 

Remarks  : 

(a)  The  inequality  (Z.4)  means  that  the  order  of  the  poles  of  A^(x) 
can  grow,  at  worst,  linearly  with  r.  This  restraining  inequality  is  the 
stronger  the  larger  k is.  The  least  upper  bound  of  the  possible  values 
for  the  number  k will  therefore  be  called  the  restraint  index  of  the  series. 

(b)  The  constant  0^  will  be  held  fixed  throughout  the  arguments 
to  come.  It  can  be  made  equal  to  zero  by  a rotation  of  the  x-plane. 
Therefore,  the  condition 

(2.7)  0Q  = 0 

does  not  restrict  the  generality  of  the  paper.  It  will  be  imposed  from 
now  on.  If  it  is  desirable  to  specify  the  other  constants,  the  corresponding 
subfamily  of  d will  be  designated  by  d(tQ,  xQ,  eQ,  6Q,  «).  Decreasing 


r 
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any  ot  these  constants  narrows  the  family,  for  simplicity,  it  will  be 
assumed  that  x^  < 1,  although  this  is  not  essential.  Then  one  must 
also  have  t^  < xn  < 1 . 

(c)  In  subdomains  of  K ^ where  |x|  > x^,  x^  independent  of  e, 

N + l 

the  error  in  2.t>)  is  uniformly  0(c  ),  as  for  ordinary  asymptotic 

series.  At  the  other  extreme,  the  best  that  can  be  said  uniformly  in  the 

k(N+1) 

whole  domain  f is  that  the  error  is  0(t„  ),  as  t_  -*  0.  In  the 

e 0 ’ 0 

I I — 2 K ^ 

subdomain  of  f-  defined  by  |x|  > tQ  e , with  «'  < k,  the  right 

(1— K' / is  ) (N+l), 

member  of  (2.  6)  becomes  0(e  ),  as  e -*  0+. 

— 

(d)  If  A(x,  e)  is  in  (7,  the  matrix  x A(x,  e),  with  k > is 

h — 

not,  unless  ( x)  = 0.  However,  e x A(x,  e),  with  h > 0,  belongs 

to  a with  a restraint  index  k1  > min(K,h/k).  This  is  most  easily 

verified  by  a geometric  interpretation  of  (2.  4):  The  points  (0,  0),  (1,  it  ), 

(2,0-  ),  etc.  in  the  (r,  <r)-plane  lie  on  or  below  the  line  r - k<t  = 0 

h -k 

and  in  the  right  half-plane.  Multiplication  by  e x translates  this 
point  set  by  the  vector  (h,  k). 

The  lemma  below  is  almost  obvious,  and  the  proof  will  be  omitted. 


Lemma  2. 1.  H A and  B are  in  <7,  then  A + B and  AB  are  in 

t .x  . e .6  .«)  with  values  of  the  constants  that  are  the  smaller  one 
' 0’  0’  0’  0’  

nl  the  corresponding  values  for  A and  B,  respectively. 

Lemma  2.  2.  If  A t tf(t,  x , c , 6Q,  k)  and  AQ(0)  is  invertible,  .then 
A_1'«  <7(t  ,x  ,cQ,6  , k)  with  constants  tQ,  xQ,  eQ  that  do  not  exceed 


V V V 
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Proof.  By  continuity,  AQ(x)  is  invertible  for  small  x.  If  t is 

small  enough,  the  relation  (2.6)  for  N = 0 implies  then  that  A(x,  e) 
is  invertible  (see  the  Remark  (c),  above).  To  calculate  an  asymptotic 
series  for  A (x,  e),  one  first  determines  recursively  a sequence  of 
matrices  B (x),  r = 0,1,  ...,  from  the  formal  relation 

00  00 

( 2,  a (X)er)(  2,  B (X)er)  = I , 
r = 0 r = 0 

whose  expansion  leads  to  the  relations 


-1 


■1  V 


B = A ; B = A ‘ > A B.,  r = 1,  2,  . . . . 
0 0 r 0 •->  r-i 

] = 0 


They  imply,  by  means  of  a simple  induction,  that 

-°rv 

(2.8)  B^(x)  = x B^(x),  B^(x)  holomorphic  in  |x|  < xQ 

Now  define  the  remainders  R (x,  e),  S^(x,  e)  by 


'N 


(x,  e)  = A(x,  e)  - X Ar(x)er 


N 

l 

r=  0 
N 


S (x,  e)  = A (x,  e)  - £ b (x)ef 

r = 0 r 

Then  the  identity  AA  * = I can  be  written  as 


N N N 

ASN  ’ 1 - <1  V “/>  - Rn  l »/  , 

o o r o r 


or,  in  consequence  of  the  definition  of  the  B , as 

r 


(2.9) 


Sm  = A'1  'l  A B e 
N . XI  r s 

r + s > N 

r,  s < N 


x , N 

- A'V,  V B er 

N-  r 
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The  summation  in  the  right  member  satisfies  the  inequality 


A (x)B  (x)er 

— j f Q 

r + s >N 
r,  s < N 


^ /-r /k~s/k  r+s  , , j — !/». 

< ! x I e < K ( I x I 


, N +1 


with  some  constant  K , by  virtue  of  (Z.  3),  (Z.  4)  and  (2.  H) . Now, 

-1  N, 

A tx,  e)  is  bounded  in  , for  0 < e < e„,  and  so  is  ).  B (x)er. 

e - 0’  r 

i-1/k  N+l  r_0 

Also,  Rjj(x,  e)  = 0 ((  |x  I e)  ) in  that  domain.  Hence,  (Z.9) 

establishes  the  same  order  of  magnitude  for  S^(x,  e),  which  proves 

the  Lemma. 

Lemma  Z.  3.  If  A(x,  e)  « t7(tQ,  xQ,  eQ,  &Q,  K)  and  t , xQ,  6Q  are  any 

positive  constants  less  than  tn,  xn,  6 then  dA(x,  e)/dx  c <7(t , xn,  e , 6. , —77; 

UUU  UUUUk+1 

and  xdA(x,  e)/dx  t £7(t  , x , e , 6 , k)  . The  asymptotic  series  for  the 

vJ  u u u 

derivative  is  the  termwise  derivative  of  the  series  for  A(x,  e).  (If 
x = 00,  one  may  replace  k/(k  + 1)  with  00. 

Proof-  One  has 


N -(0  +1)  „ 

(Z.  10)  dA(x,e)/dx  = dAQ(x)/dx  + ^ x r A (x)ef  + dR^(x,  e)/dx  , 

r=  1 


with 

A V ^ 

A (x)  = -a  A (x)  + xdA  (x)/dx  . 
r r r r 

(The  summation  in  (Z.10)  is  absent  for  N = 0.  ) Let  be  the  subdomain 

of  *!  defined  by 
e 

*e  = ^ - xo’  ^arg  ± V ’ 

where 
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- y). 


V1 


Y>'\  60  - 


f.  - 

0 


arc  sin  y 


and  y is  an  arbitrarily  small  positive  constant.  Around  a given  x t 

as  center  describe  a closed  disk  of  radius  I x I y . Then,  for  every 

t t T one  has 
x 

tj/e"  = tQ1  e*(1  " v)  - lx  ii1  “ Y)  < U I < lx|(l  + y)  < xQ(l  + y)  <xQ 


and 


I arg  | - arg  x I < arc  sin  y , 

so  that  t for  x t fl  . Hence,  by  the  Cauchy  inequality  and  (2.6), 

iye,ai/ixiv<cNixr,N+*,/K-y,tl . 

The  constant  depends  on  y. 

From  (2.  4)  it  follows  that  a + 1 < r/«  + 1 < r/-^  < r/~~,  for  r > 0. 

r K+r  k+1 

This  shows  that  the  termwise  differentiated  series  for  A(x,  e)  satisfies 
(2.4)  with  k replaced  by  k/(k  + 1).  Also,  (2.11)  implies  that  the 
analog  of  (2.6)  for  dA(x,  e)/dx  is  true  with  the  same  substitution  for  k. 
Thus,  the  statement  of  the  lemma  for  dA/dx  itself  is  proved.  The 
statement  for  xdA(x,  e)/dx  is  proved  by  going  through  the  same  steps 
after  multiplication  with  x. 

Lemma  2.4.  Let 

-a  v 

A^x)  = x r Ar(x),  r > 0 , 
where  the  ^ are  integers  such  that 

ff  < t/k 

r _ 


(2.11)  ,ldR  (x,  e)/dx|  < max 

4‘  T 
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?>  i 


lor  some  restraint  index  « iri  0 < k < *>  and  that  the  A^(x)  are 
holomorphtc  in  !x  I < x^  <1,  with  A^O)  t 0,  unless  A^(x)  = 0.  Then 

"X 

there  exists  in  every  family  <7(tQ,  XQ,  6Q,  eQ,  k)  with  tQ  < 1 a junction 

oc 

A(x.  e)  with  the  asymptotic  expansion  ^ A (x)e  . 

, r 
r=  0 

Proof.  This  lemma  is  a generalization  of  the  Borel-Ritt  Theorem  on  the 

existence  of  asymptotic  series  (see  [19]  or  Thm.  9.6  in  [IB]).  Let 

(2.12)  I A (x) I < a for  I x I < x , 

r — r — 0 


and  define  c»r(x,e)  by 

(2.13)  ar(x,e)  = 1 - exp{-a  \x  1//kc)  , 

whenever  a^  * 0.  The  positive  constant  p is  to  be  chosen  so  small 
that  the  exponent  in  braces  is  in  the  left  half  plane  for  larg  xl  < &Q. 
Then 

(2.14)  U(x,e)  I < ar\  lx|  1'^Ke)  ^ , 

for  iarg  x i < 6Q.  In  the  part  of  that  sector  where 
(2.  15)  t~lt"  < ix|  < xQ 


one  has 

° “PA 

|a  (x)or  (x,  e)e  I < ( lx  I ) e 

Therefore,  there  is  a function  A(x,  e) 


r-p  <dxrVr‘w/V'p  <t<,r'WA. 

defined  by  the  convergent  series 


00 


(2.16) 


A(x,  e)  = Yj  A (x)a  (x,  e)er  . 
r=  0 
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(Observe,  again,  that 


A(x,  e)  is  holomorphic  in  ip  for  0 < e < e 

e - O' 

one  must  choose  so  small  that  t ^e*  < x < 1) 

0 ^ 0 0 0 '* 

To  prove  that  A(x,  e)  ~ ^ A (x)e  , in  the  sense  of  the  definition 

r=0 

(2.6),  one  proceeds  as  follows. 


N N 

lA(x,  e)  £ A (x)q-  (x,  e)er|  < £ a lx  I r lexp{-a"1(x'1/,fe)_fi  | }er 
r=  0 r=0  r r 


(2. 17) 


r=  N f 1 
For  abbreviation,  set 

Then 

(2.18) 


~ “O 

+ £ a lx  I r\a  (x,  e)  |er  =:  2,  + Z . 


1 2 


Y “1/k 

f,  = x e . 


<o1A£<UI<tlA, 


for  x in  t^e*  < |x|  < x , and  also 


lx  I rer  < I ; |r  . 


Hence, 


(2.19)  Z < lilN+l  £ a exp{-a:1Re(t'h}ldr'N'1<  lilN+1C,„, 

r=0  IN 


for  |t,l  in  the  interval  (2. 18).  The  constant  C is  independent  of  t,. 
Finally, 


(2. 20) 


L2;2=N+iidN+1tkiNtl  i ur"-1-" 

r=  N+2 


iN+1 


The  proof  of  the  lemma  is  completed  by  inserting  (2.19)  and  (2.20)  into  (2.17). 
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§ 3.  The  Program. 


The  more  general  types  of  differential  equations  to  which  one  is 
led  in  the  reduction  process  to  be  performed  on  the  original  equation  (1.1) 
have  the  form 

, . h -k  dy  > 

(3.1)  C X dx  = A(X’  e)y  ’ 

where  h and  k are  integers,  h > 0,  and  A(x,  e)  is  of  class  O ■ 

For  equation  (1.  1), 

(3.  2)  h > 0,  k = 0,  K = °°  . 

Incidentally,  with  the  substitutions  x = te  and 

„ -<T  r-CT  K 

A^(x)er  = A^(x)t  r 

in  the  series  (2.2)  for  A(x,  e),  the  differential  equation  (3.1)  can  be 
interpreted  as  involving  two  independent  variables  a ' fast'  one,  t,  and 
a ' slow''  one,  x.  This  is  an  approach  used  in  much  of  the  literature, 
e.g. , in  [ 19  ] . 

The  sequence  of  transformations  through  which  one  reaches, 
eventually,  a form  amenable  to  elementary  solution  may  be  long.  To 
avoid  an  awkward  proliferation  of  the  notation,  the  previous  notation  will 
be  re-adopted  after  each  transformation.  The  simplification  of  the  problem 
reached  at  each  stage  appears  then  as  a property  that  may  be  assumed 
without  loss  of  generality. 

It  will  be  helpful  to  begin  with  a list  of  the  types  of  transformations 
to  be  performed,  to  number  them  and  to  give  them  descriptive  names. 
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Type  I (Transformation  by  power  series):  Set 

y = P(x, e)y  , 

with  P t and  P(0,e)  invertible  for  small  e > 0. 
Type  II  (Parameter  shearing): 


<*1  a a 

y = diag(e  , e , . . . , e n)y 
where  a ,a^,...ta  are  certain  rational  numbers . 
Type  III  (Independent  variable  shearing): 

P1  P2  ~ 
Y = diag(x  , x , . . . , x )y 

with  certain  rational  numbers  p^,  p^,  . . . , p . 

Type  IV  (Re-adjustment  of  the  parameter)  : 


~m 


e = e , m a positive  integer  . 

Type  V (Re-adjustment  of  the  independent  variable): 

~m 

x = x , m a positive  integer. 

Type  VI  (Eigenvalue  shifting): 

y = exp(e  ^ f \(x)xkdx}y  , 

where  \(x)  is  a holomorphic  scalar  function  (usually  a constant). 
A transformation  of  type  I changes  (3.1)  into 


(3.  3) 


h -k  dy  ,D-1  h -k  -1  dP.~ 

: x ~^=(PAP-ex  P — )y  . 
dx  dx 


Lemma  3. 1.  If  A and  P are  in  a with  restraint  indices  Kp,  if 
P(0,  e)  is  invertible  for  small  positive  e,  and  if  h > 0,  then 


-16- 


is  in  rj  with  a restraint  index  « that  satisfies 


d-!ad  h-k-ldP 
P AP  - e x P ~ 
dx 


the  inequalities 


( 3.  4)  and 


k > min(KA,Kp,  k M) 


k > min(KA,  t 1),  k) 


when  k > 0.  It  - 00 , the  middle  term  in  both  inequalities  may  be 

omitted.  If  k < 0 the  last  term  can  be  omitted  in  (3.4). 

d f h -k  dP  h -k-1  dP.  . , .. 

Proof.  Write  ex  — = e x (x~).  then  the  lemma  follows 
dx  dx 

immediately  from  Lemmas  2.1,  2.2,  2.3  and  Remark  (d)  after  Definition  2.1. 
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§4.  An  Approximate  Block-diagonalization . 


In  [ 16 1 Sibuya  proved  a fundamental  block-diagonalization  theorem 
for  differential  equations  of  the  form  (1.1).  In  this  and  the  following 
section  Sibuya's  result  is  being  extended  to  the  more  general  differential 
equation  (3.1). 


The  block  diagonalization  is  indicated  only  when 
(4. 1)  h > 0 

and  if  AQ(0)  in  (2.2)  has  at  least  two  distinct  eigenvalues.  What  follows 
applies  only  in  this  case. 

It  is  clear  that  a suitable  transformation  of  Type  I with  constant  coef- 
ficients will  change  AQ (0 ) into  its  Jordan  canonical  form.  This  transforma- 
tion does  not  affect  the  restraint  index.  According  to  the  program 
explained  in  section  3 it  may  be  stipulated  that  AQ(0)  has  already  the 
property  below. 

Property  1:  AQ(0)  is  in  Jordan  canonical  form. 

Since  Ag(0)  has  more  than  one  eigenvalue  it  can  be  written  - not 
uniquely  - in  the  partitioned  form 


(4.2) 


Aq(0) 


A>)  0 

- 0 Aq2(0)_ 


11  22 

with  AQ(0)  Aq  (0)  being  Jordan  matrices  without  a common  eigenvalue. 

There  is  an  often  proved  theorem  to  the  effect  that  there  exists  a 
holomorphic  nonsingular  matrix  P(x)  in  |x|  < xQ  (0  < xQ  < xQ)  such 
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that  P~\x)A  (x)P(x)  is  block-diagonal  in  the  same  partitioning  as 
that  in  (4.2).  (See,  e.g.,  [4]).  A transformation  of  Type  I with  this 
matrix  P replaces  AQ(x)  by  a block  diagonal  matrix,  as  described. 
Thus,  one  can  stipulate 
Property  Z: 


A0(x)  = 


"a11/  \ 

Aq(x) 

L o 


o 


where  A^  (0),  j = 1,  2,  are  in  Iordan  form  and  have  distinct  eigenvalues. 
Note  that  this  transformation  may  have  lowered  the  restraint  index 


to  h/k,  if  k > 0.  (See  Lemma  3. 1). 

Now,  a matrix  P(x,e)  < d will  be  constructed  such  that  the  new 

coefficient  matrix 


(4.  3) 


-1 . n h -k  -1  dP 
P AP  - e x P — 
dx 


of  the  differential  equation 

h -k  dz  . 

(4. 4)  e x dx  = B^X’ 6 Z 

resulting  from  (3.1)  by  the  transformation 
(4.  S)  Y = p(x»  £)z 

has  an  expansion 

00 

B(x,  e)  ~ ^ B (x)er 
r=  0 r 


which  has  all  its  coefficients  B^x)  block  diagonal  in  the  partition 
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defined  by  (4.2).  To  that  end,  one  inserts  the  series 


(4.7) 


oo 

V 


P(x,  e)~  2j  P(x)e  , 
r = 0 


as  well  as  the  series  in  (4.6)  into  (4.  3)  re-written  in  the  form 
(4.8) 


ehx~k  = AP  - PB 
dx 


and  identifies  like  powers  of  e in  the  two  members. 

The  first  of  the  recursive  sequence  of  conditions  so  obtained  is 


It  is  satisfied  by  setting 
(4.9) 


Vo  * Vo  = ° 


P0  = Bo  = Ao  • 


The  subsequent  conditions  then  become 


(4.10) 
with 

(4.11) 


A_P  - P A = B + H , r > 0 
Or  r 0 r r 


H = >.  (PB-AP)-A-x 

r LJ  u v i/  ii  r 

v + ^ = r 


-k  dPr-k 
dx 


p,  v <r 

The  last  term  in  (4.11)  is  absent  for  r < h. 

The  successive  calculation  of  the  matrices  B^,  r = 1,  2,  . . .,  from 
(4.10),  (4.11'  so  that  they  are  all  block  diagonal  follows  almost  exactly 
Sibuya's  procedure  in  (16]  as  presented  in  [18],  §§11  and  26.  One  desires 
to  have 


(4.12) 


B = 
r 


B1  1 0 

r 


b22j 

r 


, r > 0 . 
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This  can,  indeed,  be  achieved,  even  by  means  of  matrices  P of 

r 

the  special  form 


(-1.13) 


P 

r 


0 


L r 


r 


0 


as  can  be  seen  by  inserting  (4.12)  and  (4.13)  into  (4.10).  One  finds  the 


four  conditions  (cf.  [ 18  | , formula  (11 . IS)) 

,1  1 


(4.14) 


0 = B'  + H' 
r 


a>:2-«2  ■«« 

0 = B^  + 

r r 


If  B P have  already  been  calculated  for  j < r,  the  matrix  H is 

J r 

1 1 2.2. 

known.  Since  AQ  , AQ  have  no  common  eigenvalues  for  small  x, 
the  matrices  P^  and  P^  can  be  calculated  uniquely  from  the  off- 
diagonal  relations  in  (4.14).  The  matrices  b'  \ B^  are  given  directly 
by  the  diagonal  equations  in  (4.14). 

The  matrices  Pf,  so  determined  have  poles  at  x = 0 and  it 
must  still  be  shown  that  the  orders  of  these  poles  grow  at  most  linearly 


with  r,  i.e.,  that  (2. 4)  is  true  for  some  restraint  index  k. 

It  is  clear  from  (4.14)  that  if  is  already  known  to  have  a pole 

of  order  a , then  P and  B have  poles  of  at  most  order  a . By 
1 r r r 1 

induction,  it  follows  from  (4.11)  and  (4.14)  that  the  inequality  (2.4)  is 
true,  with  the  restraint  index  K of  A,  up  to  r = h.  Beyond  that  one 
has  the 

Lemma  4.1.  The  series  (4.6)  and  (4.7)  constructed  by  solving  (4.9).  (4.14) 
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satisly  ( 2 . 4 ) with  a restraint  index  k = x p which  satisfies  the  inequality 


(4.15) 


Kp  > min(xA,  7-77),  when  k>0 


A’  k + 1 


and 


Kp  > k A,  when  k < 0 , 


where  x . is  the  restraint  index  of  A. 

A 

Proof.  Assume  k > 0.  The  assertion  is  already  known  to  be  true  for 
r < h.  Assume  it  to  be  true  for  all  subscripts  below  a certain  r.  Then 
H in  (4. 11)  can  have  a pole  whose  order  is  at  most 


(4.16) 


max 
v + p,  = r 
V,  p<r 


-C. 

LK 


X-  XA 


— 1 
K * J ’ 


where  x1  is  the  right  member  of  (4. 15) . If  x1  = ~ ~ , then  x1  < xA 

and  k + 1 + ~ ~ = ^(k  + 1)  = — . Hence  the  number  in  (4.16)  equals 

x h x 

r/x 1 . If  x 1 = x . , one  has  x 1 < “TT » i • e . k+l<~,  so  that 

A — k + 1 ~ x 1 

k + 1 + ~ ~ < ~ . Again  (4. 19)  has  the  value  ~ . As  P and  B do 
x x x r r 

not  have  poles  of  higher  order  than  has  H the  lemma  is  proved  for 
k > 0.  For  k < 0 one  has,  always,  x1  = x^. 

A reference  to  Lemma  2.4  now  guarantees  the  existence  of  matrix 
functions  P(x,  e),  B(x,e)in  d,  with  restraint  index  xp  as  in  (4.15), 
such  that  P (x)  and  B^(x)  are  the  functions  in  (4.12)  and  (4.13), 
respectively.  It  does  not  follow,  however,  that  these  matrices  P and  B 
satisfy  (4.  3)  exactly  . It  is  true  that  the  two  members  of  (4.  3)  will  have 
the  same  asymptotic  expansion,  but  they  may  still  differ  by  a function 
about  which  nothing  is  known  except  that  its  asymptotic  series  has  all 
its  terms  zero. 
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According  to  the  program  of  this  paper, 
be  formulated  in  the  original  notation  as 
Property  3. 


(4. 17) 


A (x) 
r 


A*l(x)  0 

22  ’ 

L 0 A^  (x)J 


the  result  of  this  section  will 


r = 0,1,. 
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j s • Complete  Blockdiagonalization. 


Large  parts  of  this  section  resemble  the  material  in  §§14  and  26 
c>!  ! 16  j,  which  may  be  consulted  for  a more  detailed  description  of  the 
method.  See  also  [ 10  j and  [ 6 J . 

thanks  to  Lemma  2.4  and  Property  i,  above,  there  exists  a truly 
blockdiagonal  matrix  A(x,e)  in  such  that 


(5.1) 


A)j(x,  e) 


j = 1,  2 


in  the  sense  of  the  definition  in  (2.  6).  Property  3 does  not  imply  that 
A(x,  e)  is  strictly  blockdiagonal  but  if 


(5.2) 

then 

(5.3) 

with 


A(x, e)  := 


AU(x,  e) 


l 0 


~ 22, 

A (x,  e. 


A(x, e)  - A(x, e)  = : E(x, e)  , 


(5.4)  E(x,  e)  ~ 0 . 

The  purpose  of  this  section  is  to  transform  the  differential 
equation  (3.1)  - whose  coefficient  matrix  may  now  be  assumed  to  satisfy 
(5.  3)  and  (S.  4)  - into  one  with  the  strictly  blockdiagonal  matrix  A(x,  e). 
This  is  to  be  done  by  a trans formation  of  Type  I.  In  other  words,  a 
matrix  P * (J  must  be  found  such  that  the  coefficient  matrix  in  (3.  3)  is 
equal  to  A or,  equivalently,  that 
(5.  5)  ehx~kP'  = AP  - PA  . 


-24- 


In  view  of  (5.3),  (5.4)  it  is  obvious  that  (5.5)  can  be  formally  satis fied 
by  a series  o f the  form  given  in  (4.7)  with  PQ(x)  = I,  P (x)  = 0,  r > 0. 
It  must  be  shown  that  among  the  infinitely  many  functions  in  d with 
that  trivial  asymptotic  expansion  there  is  one  that  satisfies  (5.  5) 
exactly . This  will  be  done  by  following  essentially  Sibuya's  procedure 
in  [16].  (See  also  [ 18  ] , §§12,  16). 

Insert 


(5.6) 


P = 1 + 


0 P 


into  (5.  5).  By  an  easy  elimination  there  results  then  the  differential 
equation  of  Riccati  type, 


(5.7) 


h -k  dP  .22*21  *21  11  .21  121  12*21 

e x — — -A  P -PA  + A -PAP 

H v > 


"21  "12 
for  P . An  analogous  equation  holds  for  P . The  problem  now  consists 

in  proving  that  (5.7)  possesses  a solution  which  is  asymptotic  to  zero. 

The  notation  will  be  simpler  if  the  entries  of  the  rectangular  matrix 

*21 

P are  arranged  - no  matter  in  what  order  - into  an  m-dimensional  vector 
to  be  called  w.  Equation  (5.7)  can  then  be  written  as 


(5.8) 


h -k  dw  ,,  . 

ex  = f(x,w,£) 


with  a right  member  that  has  the  following  properties: 
(i)  f is  quadratic  in  w,  in  the  sense  that 


f(x,  w,  e)  = a(x,  e)  + M(x,  e)w  + g(x,  w,  e)  , 
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where  M is  an  m X m matrix  and  the  components  of  the  vector  g 
are  quadratic  forms  in  the  components  w.,  j = 1,  2,  . . . , m of  w. 

(ii)  The  vector  function  a(x, e)  is  in  d and  is  asymptotic 
to  zero.  The  reason  for  the  last  assertion  is  that  A^x.e)  ~ 0 in 
(5.7),  by  (5.2),  (5.  3)  and  (5.  4). 

(iii)  M(x,  e)  t d ■ Let 

00 

(5.10)  M(x,  e)  - ^ M (x)er 

r^O  r 

be  the  asymptotic  expansion  of  M(x,  e).  Then  MQ(0)  is  non-singular. 

11  22 

This  reflects  the  fact  that  AQ(0)  and  AQ  (0)  have  no  common  eigen- 
values, according  to  Property  2 in  section  4.  No  generality  is  lost  by 
assuming  MQ(0)  to  be  in  Jordan  form,  since  this  requires  only  a 
preliminary  linear  transformation  of  w with  a constant  coefficient 
matrix. 

(iv)  Let 


(5.11) 


m 

Vx-e'w):=  \ , Ve>(x’  c)w„v 

O', (3  = 1 


j = 1,2,  . 


, m 


be  the  components  of  g in  (5.9).  Then  g t d.  and  g ~ 0. 

jap  jap 

Again,  the  last  mentioned  fact  is  a consequence  of  the  fact  that  the 
off-diagonal  blocks  of  A(x,  e)  are  asymptotic  to  zero. 

Following  a standard  procedure,  the  differential  equation  (5.8) 
is  now  written  in  the  form 


(5.12) 


h -k  dw  ..  ...  . . 

E X diT  = M0(0)W  + p^x’  w’ 
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with 


( 5.  13)  p(x,  w,  e)  = a(x,  e)  + (M(x,  e)  - Mq(0))w  + g(x, w,  e) 

and  then  converted  into  the  equivalent  integral  equation 

h — h k. 

(5.14)  w(x,  e)  = e | exp{[  q(x)  - q(t,)  ] M (0)e  }p(C,  wit,,  e),  e)£  dl  . 

r(x) 

Here, 

k+1 

q(x)  = x /( k.  + 1),  for  k * -1,  and  q(x)  = log  x,  for  k = -1  . 

The  symbol  T(x)  stands  for  a set  of  m paths  y (x),  1 = 1,  <2,  • • • , m, 
all  ending  at  x,  each  to  be  used  for  the  integration  of  the  corresponding 
component  of  the  integrand  vector. 

The  construction  of  regions  and  paths  suitable  for  the  asymptotic 
analysis  of  the  integral  equation  (5.13)  requires  some  geometric  preparations: 
The  mapping 

1 = q(0,  with  x = q(x) 

takes  the  domain  in  the  t,-plane  into  a domain  q(/5£)  in  the  t,-plane. 

For  k > 0,  the  latter  domain  is  again  a sector  of  an  annulus  near  the  origin, 

as  in  (2.  5).  If  k = -1,  it  is  a rectangle  in  the  left  half  plane  and  for  k < -1 

a sector  of  an  annulus  in  the  left  half-plane,  close  to  the  point  at  infinity  for 

small  £ . If  X. , j = 1,  2,  ...,m,  are  the  eigenvalues  of  MQ(0)  - none  of  which 

are  zero  - there  are  finitely  many  "exceptional"  directions  in  the  t,-plane  along 

which  Re(t,k.)  is  constant  for  at  least  one  j.  Now  construct  a rhombus  ft  in 

q(C  ) according  to  the  following  rule.  If  the  bisecting  ray  of  q(;8  ) does  not 
£ £ 

have  exceptional  direction,  then  q(x^)  and  q(t^eK)  are  two  opposite 
vertices  of  ft  and  the  interior  angles  of  R at  these  vertices  must  be  so 


-27- 


small  that  no  directions  from  these  vertices  into  R are  exceptional. 

If  the  bisecting  ray  of  q(jfl^)  !§.  exceptional,  the  construction  must  be 
modified  by  rotating  the  diagonal  of  R slightly  with  respect  to  the 

real  £,-axis,  but  keeping  the  two  opposite  vertices  on  the  boundary  of  q (£  ). 

# * * # 

Next,  choose  x , t.  and  6 so  small  that  the  image  q(fl  ) 
uuu  e 

of  the  domain 


* 

6 = (x 

e 


1 « It  ^1  I ^ 

ln  E 1 U I < larg  4 le&J 


in  the  t,-plane  lies  inside  R . The  pre-image  of  R will  be  denoted  by  R. 

(See  the  figure,  below.  The  case  k > -1  has  appearance  similar  to  k > 0, 

except  that  the  sector  opens  to  the  left. ) 

For  the  remainder  of  this  proof,  x is  to  be  an  arbitrary  point 

of  R.  If  x , x are  the  vertices  of  R and  x = q(x),  then  Re(£k.) 

* 

decreases  strictly  along  one  of  the  two  directed  segments  x^x  or  x^x  in 

the  £-plane.  This  was  the  purpose  of  the  preceding  construction.  Let 

y.(x)  be  that  segment.  The  path  y^(x)  in  the  t,-plane  is  its  pre-image. 

The  existence  of  a solution  of  the  integral  equation  (5.  14)  can  be 

established  by  the  contraction  mapping  theorem.  Let  6 be  the  Banach 

space  of  vector  valued  functions  v of  x holomorphic  in  the  closure 

of  the  rhombus  R with  the  norm 

II  v ||  = max  |v(x)  | 
x < R 


where 


m 

I v(x)  I =:  £ I v (x)  I . 

J = 1 J 


For  each  e in  0 < e < eQ  the  right 
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k > 0 


■cceptional 

ase 


1 


member  in  (5.14)  defines  an  operator  on  w i 0,  which  will  be  called  C . 

e 

To  estimate  II H II,  one  first  observes  that  the  properties  (i)  through  (iv) 
imposed  on  the  function  f earlier  in  this  section,  imply  that  p,  as 
defined  in  (5.13),  satisfies  the  inequality 

(5.15)  II  P(* , w,  e)  ||  < cNt[)N+1)/K(l  + ||  w ||)2  + II  w ||  . 


Here,  N is  an  arbitrary  non-negative  integer,  and  c.,c.,  are  constants 
independent  of  e and  tQ,  for  0 < e < e . 

The  expression  for  U^w  defined  by  the  right-hand  member  of 
(5.14)  can  be  more  easily  appraised  after  the  stretching 


(5.16)  x=  xe  h = q(x)e  h,  £,  = £e  h = q(£,)e 

which  takes  it  into 


(5.17) 


(x-4)M  (0) 

O'  w = J e p(;,  W(t,,  e),  e)d;  . 

C v v 

r(x) 


I 

/ 


The  symbol  r(x)  stands  for  m paths  y.(x),  j = 1,2,  ...,m  in  the 

v v 

^-plane,  which  are  straight  segments  ending  at  x.  The  letter  t,  in  the 

v 

integrand  is,  of  course,  the  pre-image  of  t,  under  the  transformation 
( 5. 16) . Now 


(x-Om  (o)  * * ; 

e 0 e(x-;)Ae(x-C)H 


because  M^(0)  is  assumed  in  Jordan  form,  by  Property  (iii)  of  f(x,w,e). 

Here,  A = diag(k.,  X.  . . . \ ) and  H is  nilpotent.  The  entries  of 
12  m 

V V V V 

exp[  (x  - t,)H)  are  polynomials  in  x - say, 

(x-C)H  * rn 

e = {q  f(x  - ;)}  . 
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Therefore,  the  j component  of  the  vector  7 in  (5.17)  is 


(5.18) 


(x-t,)V  m 

e 1 h qj#(x  - £)Pj(£,  w(C,  e),  e)d;  . 


By  construction,  Re[  (x  - C, ) \^  j is  strictly  negative,  for  ail  x * , 
along  the  path  of  integration  in  (5.18).  The  rhombus  9 depends  on  e, 


but  the  eigenvalues  V,  and  therefore  the  directions  of  the  paths 

v v 

Vdx),  do  not.  It  follows  that  there  is  a constant  6 > 0,  independent 


of  £ , tQ  and  xQ,  so  that 


(5.19) 


m (x-£,)\.  m 


J I N 


-6  ( x-t. 


h “ C)  I < e , j = 1,  . . . ,m  , 


for  all  x in  the  image  of  9 in  the  t,-plane,  and  for  all  £,  on  •y.(x). 

The  estimates  (5.15)  and  (5.19)  suffice  to  appraise  Z w from 

e 

(5.17).  One  finds 

00 

(5.20)  II w II  < II  p(  - , w(-  , e),  e)  II  J e °pdp 

e 0 

^6'1[cNt(0N+1)A(1  + Hwll1)  + c^^llwlll  . 

The  contraction  property  of  is  proved  in  a similar  way  with 

the  help  of  the  inequality 

( 5.  21)  ||g(x,  w,  e)  - g(x,  v,  e)  II  < C ^ t ^ [ max(  ||w  ||  > ||  v II ) ) II  w - v II ) , 

which  follows  easily  from  the  property  (iv)  of  f(x,w,  e).  Formulas  (5.21) 


and  (5.13),  together  with  (5.17),  lead  to 
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t ■>.  llj£w  - o£v  j|  < 6 *[  cfY*  + CNtQN+J^kmax(  II  w II,  II  v II)  J ||w  - v II) 

in  a way  that  is  analogous  to  how  (5.20)  was  derived.  The  last  inequality 
shows  that  contracts  the  distances  in  any  given  ball,  say  in 
^ w II  < 1,  it  tg  is  chosen  sufficiently  small.  Also,  one  sees  from 
(5.  20)  that 

h (o)  II  < v < i , 

e ~ 

for  small  tQ . The  contraction  mapping  theorem  as  stated,  e.g.,  in  [ 3 J , 
Ch.  XII,  now  guarantees  the  existence  of  a unique  holomorphic  solution 
w of  (5.12). 

from  the  result  thus  obtained  one  concludes  by  a short  calculation 

that 

II w ||  < Kjjt^(t^)  for  all  N,  as  tQ  - 0 

for  x < 0 < e < Eq,  indep.  of  e.  One  needs  only  to  replace 

Hit  w ||  in  (5.20)  by  II w II  and  to  solve  the  quadratic  inequality  for  liw||. 
This  is,  of  course,  not  as  strong  as  the  claim  that  w ~ 0 in  the  sense 
of  the  definition  (2. 6),  since  t_>lx|  *ek  for  x t K and  0 < e < e„. 

However,  if  0 < k'  < k and  x is  restricted  to  the  smaller  sector 
/!'  defined  by 

E 

x t = {x/tQ  e*  < |x  I < xQ,  larg  x | < &}  , 
then  (5.23)  does  contain  the  desired  Information.  In  fact,  lx  I > tQ ^ 


means  that 


so  that 


K ” K 

Hence,  ( 5. 1 i)  is  valid  with  t replaced  by  tQe  , 

II w ||  < )M,  Km  a constant  . 

For  any  given  integer  N > 0 set  M > (N  + 1)/(k  - k')  and  thus  write 
the  last  inequality  as 

II  II  „*  N+l 
llw  II  < K^e 

As  was  assumed  to  be  less  than  one,  xq  ^ > 1»  so  indeed, 

Hwll  <(Knx-1A')N+1  for  all  N,  if  x t #1^.  The  following  theorem  has 
now  been  proved. 

Theorem  5.1.  If  A(x,  e)  has  Properties  1,  2,  3 as  stated  in  section  4, 
there  is  a transformation  of  Type  I which  takes  the  differential  equation  (_3.1j 
into  a system  with  a blockdiagonal  coefficient  matrix 


(5.24) 


AU(x,  e) 


A (x,  e). 


, AJJ(x,  e)  £ Cl,  ) - 1,  2 


By  repeated  application  of  Theorem  5.1  the  given  system  of 
differential  equations  may  be  split  into  a set  of  uncoupled  systems  of 
lower  order,  each  of  which  has  a leading  coefficient  matrix  with  only 
one  distinct  eigenvalue  at  x = 0.  This  eigenvalue  can  be  reduced  to 
zero  by  a transformation  of  Type  VI.  Accordingly,  only  differential  equa- 
tions that  have  the  property  below  need  to  be  considered  from  here  on, 
at  least  when  (4.1)  is  true. 

Property  4 . When  h > 0,  then  A (0)  is  nilpotent  and  in  Jordan  form. 
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t>.  Simplification  by  Arnold's  Theory. 


In  [ 1|,  V.  I.  Arnold  has  developed  a similarity  theory  for  holomorphic 
matrix  functions  of  one  or  several  variables,  and  in  [ 22  j I have  shown 
that  this  theory  is  useful  for  the  asymptotic  simplification  of  differential 
equations  of  the  form  (1.1). 

Here  is  a very  brief  description  of  Arnold's  canonical  form  for  the 
matrix  Ag(x),  which  is  assumed  to  have  Property  4 stated  at  the  end 
of  section  5. 

Let 

(6.1)  m>m  >*-->m,  m.  + m_  + • • • + m = n , 

12  ~ p 1 2 p ’ 

be  the  degrees  of  the  elementary  divisors  of  AQ(0),  and  assume  - this 
implies  no  loss  of  generality  - tnat  the  corresponding  diagonal  blocks  of 
Aq(0)  are  arranged  in  this  order  of  decreasing  size.  Denote  by  Aj^(O), 
p,  v = 1,  2,  . . . , p,  the  rectangular  blocks  in  the  partition  of  AQ(0) 
corresponding  to  this  Jordan  form,  and,  more  generally,  denote  by 
the  blocks  in  the  analogous  partition  of  any  n X n matrix  M. 

Arnold's  theorem,  applied  to  the  matrix  AQ(x)  says  that  it  is 
holomorphically  similar,  in  a sufficiently  small  disk  |x|  < xQ,  to  a 
matrix  Ag(x)  in  which  A^v(x)  has  non-zero  entries  in  its  last  row  only, 
when  p < v,  and  in  its  first  column  only,  when  p > v.  A matrix  of 
this  structure  will  be  called  an  Arnold  canonical  matrix.  It  will  again  be 
assumed  that  this  reduction  of  A^(x)  has  already  been  performed  by  a 
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transformation  of  Type  I independent  of  e.  In  other  words,  Property  4 
of  section  5 will  be  strengthened  to 

Property  5.  AQ(0)  is  nilpotent.  and  AQ(x)  is  in  Arnold's  canonical  form. 

This  result  can  be  improved  by  the  method  of  [22|,  as  will  now 

be  explained.  To  that  end  it  is  useful  to  introduce  certain  constant  n x n 

matrices  r,  j = 1.2,  ...,d.  The  r.  have  only  one  non-zero  entry, 
j’  J 

and  that  entry  is  equal  to  1.  It  is  located  in  the  last  row  of  the  block 
r.^v  when  p < v and  in  the  first  column  when  ^ < v.  All  such  matrices 
are  contained  in  the  set  r . . r^.  The  lemma  below,  worded  somewhat 

differently,  can  be  found  in  [1]. 

Lemma  ( . 1 (Arnold) . Corresponding  to  every  constant  complex  nxn 
matrix  H there  exists  an  n X n matrix  X and  d scalars  p., 
j = 1,  2,  . . . , d such  that 


(6.2) 

(X  and  the 


Ao(°)x  - xa0(0)  - L p.r 


j = i 


p.  are  not  unique.) 
_J 


H . 


Property  5 and  Lemma  6.1  form  the  basis  of  a further  simplification 

of  the  coefficient  matrix  A(x,  e),  at  least  when  A^(0)  is  not  zero.  (It 

2 

Aq(0)  = 0,  the  p = n,  d = n , and  any  AQ(x)  that  vanishes  at  x = 0 
is  in  Arnold's  form.)  The  method  resembles  that  of  section  4 and  also 
that  described  in  [22]  . However,  instead  of  aiming  at  matrices  B^fx)  that 
are  block-diagonal,  one  stipulates  that 
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(6.  3) 


r > 0 


d 

Vx)  = pjr(x)rj> 

00 

for  all  terms  of  the  formal  coefficient  matrix  ^ B (x)er  obtained 

r = 0 r 

from  (3.1)  by  a formal  transformation  of  Type  I.  Instead  of  (4.10)  one 
gets  now  the  recursive  conditions 

d 

(6.4)  AQ(x)Pr(x)  - Pr(x)AQ(x)  - ^ Pjr(x)rj  = Hr(x)>  r > 0 . 

Now  assume  that  H^(x)  is  already  known  to  have,  at  worst,  a 
pole  of  order  at  x = 0.  This  is  true  for  r = 1,  when  = -A^, 
by  (4.11).  Setting 

-a  -a 

(6.5)  Pf(x)  = x rPr(x),  p.f(x)  = x rp.f(x) 

-a 

in  (6.4)  and  cancelling  the  common  factor  x , one  obtains  an  equation 

v v 

for  P (x)  and  p.  (x)  which  does  have  a solution  at  x = 0,  by  Lemma  6.1. 
r jr 

By  a form  of  the  implicit  function  theorem  for  holomorphic  functions  (or 
else,  by  a simple  matrix  argument)  one  concludes  that  the  equation  for 

v v 

Pr(x),  p.^x)  be  satisfied  by  functions  holomorphic  in  a neighborhood 

of  x = 0 that  does  not  depend  on  r.  It  follows  that  (6.4)  has  solutions 
P (x),  p.  (x)  with  poles  of  order  a , at  most,  at  x = 0. 

The  induction  with  respect  to  r now  proceeds  exactly  as  in  section  4. 
The  orders  of  the  poles  of  Pf(x)  and  pjr(x)  again  satisfy  the  inequality 
(2.4)  with  a restraint  index  « that  satisfies  the  inequality  (4.15). 
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As  before,  the  result  of  this  section  will  be  stated  as  a non- 


restrictive  assumption: 

Property  6.  A^( 0)  is  nilpotent.  A^(x)  is  in  Arnold's  canonical  form,  and 

d 

(6.  6)  A (x)  = V,  p (x)r  r > 0 . 

j = l 1 1 

Property  6 prescribes  a special  form  for  all  the  coefficient  matrices 
in  the  expansion  of  A(x,  e),  unless  AQ(0)  = 0.  Lemma  2.4  then 
permits  the  construction  of  functions  p,(x,  e)  in  a with  the  corresponding 
asymptotic  expansion 


(6.7) 


j = 1,  2,  . . . , d, 


Pjolot  = 0 ■ 


The  coefficient  matrix  A(x,  e)  has  therefore  the  form 

d 

(6.8)  A(x,  e)  = Aq(0)  + Yj  p.(x,  e)rj  + E(x,  e)  , 

where 

E(x,  e)  ~ 0 

in  the  sense  of  Lemma  2.4,  i.e., 

(6.9)  |E(x,  e)  I < CN(  lx  | 1//Ke)N+1 

for  all  N and  for  x t £ , 0 < e < e . 

e — o 

On  the  basis  of  the  theory  developed,  so  far,  it  cannot  be  asserted 
that  E(x,  e)  may  be  taken  as  identically  zero,  as  was  proved  to  be 
possible  in  section  5 for  the  formal  simplification  of  section  4. 
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§ 7 . Shearing  and  Re-adiustment. 


Following  a technique  used  by  Turrittin  [ 17  J and,  in  a form  more 
similar  to  the  one  here,  in  Iwano  and  Sibuya  ( 10 ) , Iwano  [6],  [ 7 ],  a 
further  reduction  of  the  problem  will  now  be  achieved  by  means  of  a 
combination  of  shearing  transformations  of  Types  II  and  III  (see  § 3). 
This  is  necessary  and  meaningful  only  if 
(7.1)  n > 1 , 

since  a scalar  first  order  linear  differential  equation  can  be  trivially 
solved.  First,  a few  convenient  symbols  will  be  introduced: 

Let 


(7.2)  p (x)=  i p xff,  |x|<xQ 

<x  = 0 

(cf.  (6.5)).  Denote  by  P the  set  of  ordered  triples  (j,r,  v)  for  which 

p * 0.  Next,  call  (u.  , v ) the  row  and  column  numbers,  respectively, 
jro  j J 

of  the  one  non-zero  entry  of  the  matrix  r . Finally,  define  12(a)  by 

(7.  3)  Q(a)  :=  diag(l,  a,  a2,  ... , a°  *)  . 

The  shearing  transformation  (Types  II  and  III  combined) 

. a (3 

(7.  4)  y = P(e  x )z  , 

with  a and  p to  be  suitably  determined,  takes  the  differential  equation 
(3.1),  in  which  A (x,  e)  has  now  the  form  (6.8),  into 


(7.  5) 


h-a  -k-p  dz 
e x — 
dx 


B(x,  e)z  , 


where 
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(7.6) 


d -6  a -6  (3 

B(x,e)  = A (0)  + V p (x,  e)e  J x 1 F 
j^l  J J 


h-o  -k-13-1  , _v. 

-e  x diag(0,  ^3,  . . . , (n  - 1)(3)  + L (x,  c)  . 


Here, 

(7.7)  6 = u - v . + 1 

) v J 

and  E (x,  e)  ~ 0.  If  the  p.(x,  e)  in  (7.6)  are  replaced  by  their  series 
expansions  from  (6.5)  and  (7.2)  the  terms  of  the  resulting  summation  are 

r-6  .a  a -a  -6,(3 

(7.8)  p e j x f J , (j,  r,  <0  t P . 

jra 

The  triples  of  P with  6^  > 0 correspond  to  terms  of  (7.6)  that  occur 
on  or  below  the  main  diagonal,  while  all  terms  with  6 < 0 are  above 
the  main  diagonal.  The  subset  of  P with  6^  > 0 will  be  called 
P ("L"  for  "lower").  Note  that,  since  6 ,n(0)  - 0,  no  triple  ( j , 0 , 0 ) 
is  in  P. 


The  parameters  a,  p are  to  be  chosen  so  that  one  term  in  (7.8) 
becomes  a constant.  For  a convenient  description  of  the  proper  choices 
for  a and  (3  the  pointsets 

(7.9)  u = r/6  . , v = (o  - <rr)/6^,  ( j , r,  a)  i P 

in  an  auxiliary  (u,  v)-plane  is  introduced.  These  points  do  not  exist 
for  6.  = 0,  unless  the  (u,  v)-plane  is  interpreted  as  a projective  plane 
and  those  points  are  taken  at  infinity  in  the  direction  indicated  by  (7.9). 
For  the  construction  that  follows  these  points  are  irrelevant,  in  any  case 
The  set  defined  by  (7.9)  will  be  called  0,  and  is  the  subset 


-39- 


corresponding  to  6 > 1.  Because  of  the  inequality  (Z.  4),  the  set  Oj 
lies  in  the  closure  of  the  sector  S bounded  below  by  the  line 
(7.10)  u+kv=0 

and  to  the  left  by  the  v-axis.  Since  (j,0,0)  is  not  in  p,  the  origin 
of  the  (u,  v)-plane  is  not  in  Q.  Also,  the  fact  that  J 6^  | <n  implies 
that  neither  the  set  Q itself,  nor  its  projection  on  the  coordinate 
axes  has  accumulation  points  in  the  finite  plane.  The  set  Q - 
lies  in  the  reflection  of  the  sector  S in  the  origin. 

To  take  into  account  the  term  with  the  powers  e*1  ^x  k ^ in 
(7.6),  the  point  with  coordinates  u = h,  v = -k  - 1 is  added  to 

* 

except  when  (3-0,  and  the  so  augmented  set  will  be  called  Q^. 

If  « > h/(k  4 1)  > 0,  the  line  (7.10)  must  be  replaced  by 


u 4 


h 

k 4 1 v 


0 


'i* 

and  S must  be  accordingly  increased  to  a larger  sector  S . 

The  determination  of  a and  p that  will  now  be  described  is 

related  to  the  use  of  the  Newton-Puiseux  polygon  in  several  other  papers 

on  the  asymptotic  theory  of  differential  equation;  for  instance  in  [ 10],  [6], 

| 11  ] . Let  n be  the  boundary  polygon  of  the  convex  hull  of  Qj  and 

II  the  boundary  polygon  of  the  convex  hull  of  Q - Q . Observe  that 
U L. 

* # 

II L is  never  empty,  since  always  contains  the  point  (h, -k  - 1). 

However  0 - may  exceptionally  be  empty. 

* 

Among  the  left-most  vertices  of  there  is  one  with  the  smallest 

ordinate.  Let  this  vertex  be  (uQ,  VQ)-  Two  ^uite  different  cases  must  now 
be  distinguished. 
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Case  1 . u0  < h. 

Take  a uQ,  (3  = vQ. 

Case  2.  uQ  = h (uQ  > h is  impossible). 

Take  a = h , (3  = 0. 

The  second  case  will  be  postponed,  for  the  time  being.  In  the 

tirst  case,  there  is  now,  indeed,  at  least  one  term  in  (7  .8)  that  has  both 

exponents  equal  to  zero,  and  for  which  (j,r,  <r)  ( Pf.  Also,  for  all  triples 

in  P one  has  r - 6 a > 0.  It  must  still  be  shown  that  the  new  non-zero 

) ~ 

exponents  of  £ and  x in  (7.8)  satisfy  an  inequality  such  as  (2.4), 
i.e.,  that  for  a suitable  positive  number  k'  one  has 
(7. 11)  cr.-(r+6.(3<(r  - 6 a)  /«'  . 

To  find  such  a k',  observe  that  the  point  set  Q: 

a - a 

(f- - — 51  - P),  (j,r,<r)  t P, 

’j  j 

is  obtained  from  the  set  Q by  the  translation  (-a,-p).  The  point  u^,  v^ 

of  Qj  becomes  the  origin,  which  is  not  counted,  since  its  contribution 

is  combined  with  AQ(0).  Therefore,  there  exist  infinitely  many  straight 

lines  through  the  origin  with  non-positive  slope  which  separate  the 

translate  Q of  0,  from  the  translate  of  Q - Q.  . This  remains  true 
L L ^ 

~ /v  3{C 

if  Qj  is  replaced  by  Q^,  the  set  obtained  by  adding  the  point 
(h  - a,  -k  - 1 - P).  Let 

u + k'v  = 0 with  0 < k'  < 00 
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one  of  those  lines.  The  number  k'  cannot  exceed  k.  Tor  later  use 
the  stronger  restriction 

(7 . 12)  K ' < K 


will  be  imposed. 


For  all  points  of  0,  one  has 

L V 

(T  - (J  . 

f- ~ ~ P)  > 0,  b > 0 , 

j j J 

which  is  equivalent  with  (7.11)  for  those  points.  For  points  in  Q - 

the  number  6^  is  negative.  As  these  points  lie  below  the  line  u + k'v  = 0 

the  inequality  (7.11)  is  again  true.  It  is  also  valid  for  (h  - o,  -k  - 1 - fl). 

The  transformed  differential  equation  (7.  S)  is  not  yet  of  the  form 
(3.1),  because  a and  p are,  in  general,  non-integral  rational  numbers 
with  least  denominator  not  exceeding  n.  By  appropriate  re-adjustment 
transformations  of  Types  V and  VI,  say, 

(7.13)  x = x , e = e , 

the  differential  equation  can  be  restored  to  a form  involving  only  integral 
powers  of  the  independent  variable  and  the  parameter.  The  new  differential 
equation  is  now 


(7.14) 

with 

(7.15) 


j(h-or)m ~-(k+p+l)q+l  !U  = B(x,  c)z  , 


B(x,  e)  = qB(xq,em)  , 
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B(x,  e)  being  as  in  (7. 6).  The  right  member  of  (7 . 6),  when  expressed 

in  terms  of  x,  e,  can  be  written  as  a power  series 


(7.16) 


B~(x)er  , 


whose  coefficients,  thanks  to  the  construction  of  the  shearing  transforma- 


tion, have  the  form 


(7.17) 


B~(x)  = X r B.(x)  , 


~ ~ | . i/Q 

where  B~(x)  is  holomorphic  in  |x|<xQ  and 


(7.18) 


a ~ < r/Z  , K = mx'/g 


It  remains  to  verify  that 


(7.19) 


B(x,  e)  ~ )]qB~(x)er 


in  the  sense  of  the  definition  in  (2.6),  i.e.,  that  the  N partial  sum 
B(N)(x,  e)  of  the  series  in  (7.19)  satisfies  an  inequality  of  the  form 


(7. 20) 


(x,e)  - B(N)(x,  e)  I < CN(  lxl"1/K?) 


l/S-.N+l 


for  x t 6~,  0 < e < eAm.  Here  is  defined  by  inequalities  of  the  form 


(7.  21) 


t”1eK  < lx  I <x^e,  Iargx|<&0/m 


To  prove  (7. 20),  denote  by  p(.M)(x,  e)  the  Mth  partial  sum  of  the 
series  for  p (x,  e)  in  (6.7).  Then  there  is  a constant  cM  such  that 


(7.22) 


<M>fx.  e)  I <c„dxrlAe)M+l 


I Pj(x,  e)  - p.  (x,  e)  I < cM( 
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for  1 < j < d,  x t ,0  < e < e . If  (x,  e)  is  the  function 


obtained  by  replacing  p.(x,  e)  in  the  right  member  of  (7. 6)  by  p ^ (x,  e), 
one  has,  therefore 

(7.2  3)  I B(x,  e)  - BtM)(x,  e)|  < cM(  |xr1Ae)M+1e"n“|xrn,i 

in  the  same  domain.  The  integer  M can  be  chosen  so  large  that  all 
terms  of  the  partial  sum  B^  (x,  e)  occur  in  qB^ ^(xq,  em) . Then 

(7.2-1)  S(x,e)  - B(N)(x,E)  = q(B(xq,Em)  - B(M)(xq,£m))  + Q(M)(£,e)  , 

where  (x,  e)  has  the  form 


(7.25) 


Q(M*(x,e)  = q ~(x)x  rer 


Mm 


r = N+l 


v l / 

with  q~(x)  holomorphic  in  lx)  < xQ7  . From  (7.23),  (7.24),  (7.25) 
one  gets,  using  also  (7.18),  that 


(7.26) 


I B(x,  e)  - B^N\x,  e)  | < qcM(  |x  | 1//kc)M+1e  na|x|  n(i 

tk..(i;r1/::-'N+1 


M' 


e ) , (kw  a constant)  . 

M 


This  is  valid  for  x t ~ , 0 < e < provided  tQ  < tg/q.  Now, 

(7.26)  can  be  written 

(7.27)  |B(x,e)  - B(N)(x,e)  I < ( Ixr1//K  e)N+1l  kM  + qcM  Ixl  Me  M 


J/q 


with 


A__  = l*  + D)  + n(3 Jq  - k \n  + 1) 
M 

B, , = [(M  + 1)  - na  | m - (N  + 1)  . 

M 
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By  (7.21), 

l~  | ~-1~k  , ^M~^M  ~(M+l)m(l-K’/t<)-nqt<  [3-nma 

ix  i e < (tQ  e ) e = tQ  e' 

In  the  last  equality  the  second  formula  of  (7.18)  has  also  been  used. 

Thanks  to  the  inequality  (7.12),  M can  be  chosen  so  large  that 

the  exponent  of  e in  the  last  formula  is  positive.  With  M so  chosen 

^ 1 /m 

the  expression  in  brackets  in  (7.27)  is  bounded  for  x t £~,  0 < e < 

This  completes  the  proof  of  (7.  20),  i.e.,  of  (7.19).  With  the  definition 

D(x)  = x 'k  f,H1^qdiag(0,  (3,  . . . , (n  - l)p) 
and  a return  to  the  notation  of  (3.1),  i.e.,  replacing  x,  e,  h - a,  k - |3, 
B(x,  e)  by 

(7. 28)  x,  e,  h,  k,  A(x, e)  , 

respectively,  the  result  of  this  section  can  be  stated  in  the  form  of  a 
non-restrictive  assumption,  as  follows. 

Property  7.  H h > 0,  the  coefficient  matrix  A (x,  e)  in  1 3. 1)  has  the  foi 

d h 

(7.29)  A(x,  e)  = AQ(0)  + £ p.(x,  e)r.  + D(x)e  + E(x,  e)  . 


(7.  30) 


V0)  * A01  + *02 


with  A nilpotent  and  in  Iordan  form,  and  AQ^  is  a lower  triangular 


matrix  of  the  form 


(7.  31) 


_ S’ 


aojrj 


4S- 


The  I'  are  the  matrices  defined  earlier  in  this  section.  The  r>  (x.c) 
j J 

are  in  class  c7,  with  lim  p.(0,e)  = 0,  and  E(x,  e)  ~ 0.  The  matrix 

e -*0  + J 

D(x)  is  diagonal. 

In  general,  the  constant  lead  matrix  Aq(0)  in  (7.29)  has  more 
than  one  eigenvalue.  The  reduction  process  can  then  be  continued  as 
described  in  section  IS.  However,  the  possibility  that  Ag(0)  has  again 
only  one  eigenvalue  must  also  be  explored.  This  is  the  subject  of  the 
next  section. 

Now  Case  2 above,  in  which  the  shearing  exponents  are  a = h, 
p = 0,  remains  to  be  considered.  If  that  situation  arises  e disappears 
from  the  left  member  of  (7.5).  The  asymptotic  solution  of  such 
differential  equations  - in  the  original  notation  of  (3.  1)  - can  be  found 
in  section  10. 
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§ 8.  An  Exceptional  Case:  All  Eigenvalues  Equal  After  Shearing. 

In  this  section  it  will  be  assumed  that  the  lead  matrix  A^(0j  in 
(7.29)  has  only  one  distinct  eigenvalue.  Then  the  block-diagonalization 
process  of  sections  4 and  5 is  omitted  from  the  sequence  of  operations 
in  the  reduction  process.  The  steps  to  be  performed  next  are:  Reduction 
to  Arnold's  form,  shearing,  re-adjustment  ( from  x,  e to  x,  e), 

a (5 

multiplication  by  e x . 


(8. 1)  £hx  k dt  = A^X’  E^y 

be  the  problem  so  obtained,  and  assume  that  the  new  lead  matrix  Aq(0) 
has  again  only  one  distinct  eigenvalue.  Here  A(x,  e)  is  the  matrix 

<v 

called  B(x,  e)  in  (7.19).  It  has  the  following  structure: 


(8.2) 


d 

A(x,  e)  = qA  (0)  + ^ p (x,  e)T 
j ~ 0 J 1 


V°)  = A01  + ,^o  Kjrj 


The  matrix  AQ^  is  as  in  (7 . 30),  i.e., 


(8.4) 


A = H 6 H © •••  ©H  , 

01  1 2 p 


where  H is  a "shifting  matrix"  of  order  m (see  (6.1)).  (By  definition, 
) J 

a shifting  matrix  of  order  greater  than  one  is  one  that  has  a line  of  Is 
above  the  diagonal  and  zero  entries  elsewhere.  A shifting  matrix  of  order 


one  is  zero.)  The  * are  constants.  They  are  zero  for  all  j with  p.  < v., 

) J J 
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The  eigenvalues  of  A^( 0 ) , which  is  a lower  block-triangular 
matrix  in  the  partition  induced  by  the  blocks  in  (8.4),  are  the  eigenvalues 
of  its  diagonal  blocks.  By  assumption,  there  is  only  one  distinct  eigen- 
value, kQ,  present.  Let  T , s = 1,  2,  . . . , p,  be  constant  mg  X mg 
matrices  such  that  a transformation  of  Type  I with  the  block  diagonal 
coefficient  matrix 

(8.  s)  t = ta  ®t2  e — e>  Tp 

takes  the  diagonal  blocks  of  AQ(0)  into  their  Jordan  form 

(8. 6)  kQIs  + Hg,  s = 1,  2, . . . , p 

(I  is  the  identity  matrix  of  order  m ).  Such  a transformation  takes 
s s 

each  block  A^V(x,  p)  of  A(x,  e)  into  T^A ^ (x,e)T^,  so  that,  in 
particular,  any  block  AgV(0)  that  is  a zero  matrix  remains  a zero  matrix 
after  the  transformation.  After  that,  the  entries  in  the  main  diagonal 

can  be  removed  by  a transformation  of  Type  VI.  The  new  coefficient  matrix 
A(x,  e)  so  obtained  has  a lead  matrix  of  the  form 

(8.7)  V0)  ' *01  **02  ' 

with  A as  in  (8.4)  and  A strictly  lower  blocktriangular. 

U 1 C* 

Next,  Aq(0)  is  changed  into  its  Jordan  form  by  another  transformation 
of  Type  I with  a constant  matrix.  This  is  not  necessary  if 


(8.8)  A02 

Otherwise,  the  leadmatrix  is  now  again  a blockdiagonal  matrix  composed 


I 
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of  blocks  such  as  (8.6),  but,  the  unique  eigenvalue  will  generally  not 
be  the  same  as  in  the  analogous  transformation  of  AQ(0).  Moreover, 

and  this  is  decisive,  a Lemma  of  Turrittin  in  [ 17  J,  (proved  somewhat 
differently  in  [ 18)  as  Lemma  19. -4),  says  that  not  all  blocks  will  have  the 
same  size  as  before  in  (8.6).  In  fact  the  Lemma  assures  that  after  a 
finite  number  of  such  sequences  of  transformations,  the  lead  matrix 
will  consist  of  one  block  only. 

If  (8.8),  is  true,  (8.7)  shows  that  AQ(0)  is  itself  a nilpotent 
Jordan  matrix  and  the  next  shearing  can  be  performed  immediately.  No 
change  in  the  block  structure  occurs  then. 

Thus,  under  the  assumption  of  this  section  that  all  shearings 
introduce  lead-matrices  with  only  one  distinct  eigenvalue,  one  is  eventually 
led  to  a problem  with  the  additional  property  that  all  nonzero  entries 
introduced  by  subsequent  shearings  are  in  the  last  rows  of  the  diagonal 
blocks.  If  only  one  block  is  present,  this  is,  of  course,  no  further 
restriction. 

The  shearing,  in  this  situation,  introduces  exponents  a and  (3 
that  are,  by  necessity,  integers,  because  there  must  appear  nonzero 
entries  on  the  main  diagonal  of  the  lead  matrix  after  the  shearing,  other- 
wise the  new  unique  eigenvalue  would  have  to  be  zero,  contrary  to  the 
fact  that  at  least  one  nonzero  entry  has  been  added  by  the  shearing.  Now, 
the  shearing  introduces  nonzero  diagonal  entries  only  if  6 in  (7.7) 


-49- 


is  equal  to  one  in  the  definition  of  (uQ,  v ) - (a,  (i).  In  that  case  a,  p 
are,  indeed,  integers. 

It  follows  that  no  re-adjustment  is  necessary:  x = x,  e = e. 

Therefore  h = h - a.  If  o > 0 occurs  h times,  at  most,  in  the 
repetitions  of  the  reduction,  a problem  without  a power  of  e in  the  left 
member  has  been  reached  and  one  proceeds  as  in  section  10. 

There  remains  now  one  final  possibility:  The  shearing  exponent  a, 
in  the  situation  described  in  the  last  two  paragraphs  is  forever  zero  in 
each  repetition  of  the  reduction  process.  This  case  is  taken  care  of  in 
the  lemma  below. 

Lemma  8.  1.  Assume  that  the  given  differential  equation  has  the  form 
described  in  (8.1)  through  (8.4)  and  that  A^fO)  has  no  nonzero  entries 
outside  of  the  diagonal  blocks.  Assume  that,  in  addition,  the  sequence 
of  operations:  Reduction  to  Arnold's  form,  eigenvalue  shifting,  shearing. 
multiplication  by  the  appropriate  power  of  x,  always  leads  to  a problem 
of  the  same  type  with  the  same  value  of  h.  Then  the  matrix  AQ(x)  has. 
identically  in  x,  exactly  one  eigenvalue. 

Proof.  To  simplify  the  notation  the  tildes  will  be  omitted  from  the  notation. 

The  shearing  transformations,  as  well  as  the  reduction  to  Arnold's  form, 
subject  A(j(x)  to  similarity  transformations  and,  therefore,  do  not  change 
its  eigenvalues.  Let  A^m^(x)  be  what  becomes  of  A()(x)  after  m 
iterations  of  the  reducing  cycle,  with  A^(x)  = A^(x).  If  p^m^(x),  ..  . , (x) 

are  the  eigenvalues  of  A^^x),  then,  by  assumption, 
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18.9) 


n . 


lim  ;m,(x)  = t 0,  j = 1,  • . • , 

x - 0 J 


th 


Let  13  the  shearing  exponent  in  the  m iteration.  Then  one  has  also 
m 


(8.10) 


(mfl).  , ^m,  (m),  . . (m).  n . 

Pj  (x)  = x (p  '(x)  - K'Q  ),  m - 0,1,...  , 


i.e.  , 


(m).  , (m+1),  . m (m) 

Pj  (x)  = p^  (x)x  + 


Hence, 


Prt+P, 


(°),  X (1),  v P0  (0)  (2),  . K0  ' 1 (1)  ' 0 (0) 

ij  (x)  = p^  '(x)x  + \Q  = pj  (x)x  + K0  ‘ 


By  repetition  of  this  procedure  one  proves  that,  for  arbitrarily  large  /, 

t -1  t 

(y.  11)  ^'(x)  - x<«  ♦ K^x1’0  ♦ ...  ♦ X<'>xS  = ° “*  ♦ O(x^0"S)  . 

On  the  other  hand,  the  functions  p|°\x),  as  eigenvalues  of 
A^(x),  which  is  holomorphic  at  x = 0,  possess  convergent  expansions 
in  ascending  fractional  powers  of  x.  By  the  uniqueness  of  power  series 
expansions  this  series  must  be  the  one  constructed  in  (8.11),  i.e., 

m-1 

«,  L es 

,<°>(x)  = l V<m,xs  = 0 . 

J m = 0 

This  series  is  independent  of  j,  which  proves  the  lemma. 


In  the  situation  described  by  the  lemma  above,  i.e.  , if  A (x)  has, 
identically  in  x,  exactly  one  eigenvalue  \Q(x)  an  eigenvalue- 
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shitting  transformation  of  Typo  VI  with  k(x)  “ produces  3 

similar  problem  with  a new  function  (, x ) which  is  identically  nilpotent. 

Aq(x)  may  already  be  assumed  to  be  in  Arnold's  form.  Therefore  its 
diagonal  blocks  are  shifting  matrices.  After  the  next  shearing,  either 
the  new  lead  matrix  Ag(0)  has  nonzero  entries  in  the  blocks 
below  the  diagonal  blocks  or  the  nonzero  entries  originated  from 
matrices  A^(x)  with  r > 0,  in  which  case  a > 0.  In  either  case 
situations  have  now  been  reached  which  were  discussed  before,  and  the 
whole  process  must  terminate  after  a finite  number  of  steps. 

The  theorem  below  summarizes  what  has  been  achieved  so  far. 
Theorem  8.1.  By  a finite  number  of  transformations  of  Types  I to  VI. 
described  in  $ 3,  every  differential  equation  of  the  form  (3.1)  can  be 
reduced  to  a finite  number  of  equations  of  the  same  form  for  which  either 
n = 1 or  h = 0. 
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§ 9 . Asymptotic  Solution  of  the  Scalar  Equation. 


If  (3.1)  is  a scalar  equation,  one  of  its  solutions  is 

-h  X k 

(9.1)  y(x,  e)  = exp{e  | £ A(£,e)d£}  . 

xo 

In  this  section  the  asymptotic  nature  of  this  solution  will  be  investigated. 

When  k < 0,  the  integral  in  (9.1)  is  not  always  in  class  <7. 
Therefore  the  following  wider  classes  are  introduced. 

A' 

Definition  9.1.  A function  A(x,  e)  is  said  to  be  in  class  (7  , if  it  has 

v 

all  the  properties  of  class  (7  except  that  A{x)  need  be  holomorphic  in 

0 < x < x,  |arg  x I < 6 

only.  But  A (x)  must  be  bounded  in  that  domain.  A function  A(x,  e) 

jJ-  :ic  # 

such  that  xPA(x,e)  t <7  for  some  p>0  is  said  to  belong  to.  .class  O • 
The  asymptotic  relation  (2.2)  is  defined  in  these  wider  classes 
exactly  as  in  Definition  2.1,  and  the  lemmas  of  Section  2 remain  valid 
in  them . 

The  integrals 

x 

(9.2)  B (x)  = j £ Ar(£)d£  , 

X0 

with  A (x)  as  in  (2.2),  (2.  3)  have  convergent  ascending  series  in  powers 
of  x,  except  for  the  possible  occurrence  of  a logarithmic  term.  The 
series  may  begin  with  a negative  power  of  x. 


Lemma  9 . 1 • Let 


and 

u X . 

(9-4) 

F(x,e)  = e"  J 4 A(|,e)d&  - Q(x,e) 

. 

xo 

Then 

* 

F(x,  e)  t a , and 

(9.  5) 

00 

F(x,  e)  ~ Yj  b h • 

r=  h+1  r 

00 

Proof. 

It  is  convenient  to  re-write  the  series  in  (9.  5)  as 

Y cr(*)e 

r=  0 

(9.6) 

( 

ZQ(x)  = 0;  Cf(x)  = Bf+h(x),  r > 0 . 

Then 

(9.7) 

lcf(x)|  < { 

const.  , for  k + 1 - 

k +*-0  r 

const.  1 x | , for  k + 1 - 

o > 0 

r+h 

a < 0 

r+h 

^const.  Ilog  x|  , for  k + 1 - 

o = 0 , 

r+h 

because  of  (9 . 2)  and  (2.  3).  As  CQ(x)  = 0,  the  inequalities  (2.  4)  and 

(9. 7)  imply  the  relations 

v 

> v v "ar  " 

(9.8)  C^x)  = Cr(x)x  , C^x)  bounded 

V V V 

(9-9)  of  < t/k,  0 < k < x < » . 

If  k < h/(k  + 1)  one  has  « = k. 

To  prove  (9.  9),  recall  that  by  (2.  6), 

N 

(9.10)  A(x,e)  = ]]  Ar(x)er  + (x  1^Ke)N+iRN(x,  e)  , 

r=  0 


with 
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Let  M > N > 0, 


with  R^(x,  e)  uniformly  bounded  in  tor  0 < e < eQ. 

replace  N by  M in  (9 . 10)  and  insert  it  into  (9  • 4) . Referring  to  (9 . 2) 
and  (9.6)  one  then  obtains  the  relation 


N Mi 

F(x,  e)  = Yj  C (x)er  + Y C (x)er 
r = 0 r=  N+l 


(9.11) 


M+l-h 


f RMU,e)6k-(M+1)Ad5  . 


From  (9.6)  and  (9.9)  one  concl  udes  that  the  second  summation  in  the 
right  hand  member  of  (9.11)  satisfies  the  inequality 


(9.12) 


I Y C (x)eF  I < const.  |x  Ke 
r=  N+l  r 


-1/k  iN+1 


If  M is  so  large  that  k - (M  + 1)A  < -1,  the  last  terms  in  (9.11)  is 
less  than 


(9.13) 


, M+l-h  | ik+l-( M+1)/k 
const,  e |x  I 


in  modulus,  i.e.,  less  than 


(9.14)  const,  e 


M+l-h-(N+l).  |k+l-(M+l)A+(N+l)A  | -1A  | N+l 


One  verifies  immediately  that  for 


t^e"  < lx  I <xQ 


the  last  expression  is  0(  lx  ^ e l^+^),  provided  M is  taken  so  large 
that  the  exponent  k + 1 - (M  + 1)A  + (N  + 1)A  in  (9.14)  is  negative. 


Insertion  of  these  estimates  into  (9.11)  shows  that 


i 
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I 


F(x,  e)  = Y C (x)ef  + 0(  lx  ^Ke  I 
r=  0 r 


for  0 < e < e^,  x t £ , which  completes  the  proof  of  the  lemma. 


From  the  expansion 


00 

F(x)  ~ C (x)er,  C (x)hO, 
r-  0 

which  is  the  same  as  (9.  5)  in  simpler  notation,  one  proves,  just  as  for 

~ jjc 

ordinary  asymptotic  power  series,  that  Y(x,  e)  :=  exp(F(x,  e))  also  is  in  O ■ 
Thus,  the  theorem  below  has  been  proved. 

Theorem  9.1.  For  n = 1,  the  differential  equation  (3.1)  has  a solution 
of  the  form 

(9.15)  Y(x,  e)  = Y(x,  e)e^X,£^ 

with  Y(x,  e)  t a'  and  Q(x,  e)  a polynomial  of  degree  h in  e *,  as. 
described  in  (9.2).  (9.3). 

Remark.  It  is  easy  to  analyze  the  function  Q(x,  e)  in  more  detail.  This 
will  not  be  elaborated  upon  here.  Since  the  result  depends  decisively 
on  whether  k > -1,  = -1  or  < -1  the  description  is  a little  tedious. 

In  all  cases  the  first  term,  B0(x)e_h,  dominates  the  others  in  at 

least  for  large  tQ. 
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§10.  Asymptotic  Solution  when  h - 0. 


The  problem  is  to  study  the  solutions  ol  a system  of  the  form 


(10.1)  x ^ dx  = A(x,e)y,  A( x , e ) < a 

for  small  e in  the  region  . Although  the  character  of  the  solution 
depends  considerably  on  whether  k > -1,  k = -1  or  k < -1,  the  discussion 
will  be  unified  as  much  as  possible. 

Formally.  (10.1)  always  admits  matrix  solutions  of  the  form 

00 

\ ' r 

Ij  Y (x)e  . The  Y (x)  must  only  satisfy  the  differential  equations 
r-0  r r 

. dY 

(10.2)  x"K  — f = A (x)Y  + Yj  A (x)Y  (x),  r = 0,  1 

dx  ur  u v 

p+v  = r 

^>0 


For  r = 0,  that  equation  is 


(10.  3) 


*'k  = Vx,Yo  ■ 


The  nature  of  its  solutions  near  x = 0 is  described  by  the  next  lemma. 
Lemma  10.1.  The  differential  equation  (10.  3)  possesses  a particular 
fundamental  matrix  solution  of  the  form 


(10.4)  Yq(x)  = YQ(x)xGeQ(X^ 

with  the  following  properties- 

(i)  YQ(x)  has  an  asymptotic  series  in  powers  of  x^13,  as 

x-*0  in  I arg  x I < 6Q. 

det  YQ(x)  4 0 for  x 4 0 and  Yg(0)  ^ 0 , 
but  det  Yq(0)  may  be  zero.  The  number  p is  a positive  integer. 
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(>»)  Q(x)  is  a diagonal  matrix.  If  it  is  not  identically  zero,  it  is 
a polynomial  in  x ^ without  constant  term. 

(iii)  G is  a constant  matrix,  which  commutes  with  Q(x). 

(iv)  I_f  k = -1,  then  p = 1,  Q(x)  = 0 end  Y^(x)  is  holomorphic 
at  x = 0.  If  k > 0,  then  G = 0,  as  well,  and  det  Yq(0)  / 0. 

This  lemma  summarizes  well  known  classical  results  which  can,  e.g., 
be  found  in  [ 1 «|  Chapters  I-V. 

Two  preliminary  transformation  will  now  be  performed.  While  not 
necessary,  they  are  helpful  in  that  they  simplify  the  subsequent  analysis. 
First,  set 

(10.5)  x = x1/p 

in  (10.1).  This  produces  a differential  equation  in  which  k is  replaced 
by  p(k  + 1)  - 1 and  p becomes  equal  to  1.  As  often  before  in  this 
paper,  the  notation  will  be  changed  back  to  the  original  one.  Then  the 
non -restrictive  condition  below  can  now  be  imposed. 

Property  8.  In  Lemma  10.1.  one  has  p = 1. 

Next,  set 

(10.6)  T(x)  = Y^(x)x^  , 

for  abbreviation,  and  perform  the  transformation 
(10.  7)  y = T(x)z 

in  (10.1).  There  exist  real  numbers  g^ , g^  (they  are  not  unique)  such  that 
g j g 

(10.8)  T(x)x  , T (x)x  ^ are  bounded  in  0 < | x I t xQ,  |arg  x|  < . 
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The  differential  equation  (10.1)  is  now  changed  into 


(10.9) 


-k  dz  D , . 

x — = B (x,  e)z  . 
dx 


Formally,  the  series  in  the  expansion 


(10. 10) 


A(x,  e)  ~ Yj  Ar(x)eP,  e - 0+  , 


i » ~ r 

becomes  2j  B (x)e  , with 
r-  0 r 


(10.11) 


B (x)  = T (x)Ar(x)T(x),  r > 0 . 


The  new  leading  coefficient  BQ(x)  is  more  easily  calculated  indirectly: 

OO 

In  analogy  to  the  formal  series  solution  ^ Y (x)e  the  differential 

r=0  r 

00 

equation  (10.9)  has  a series  solution  ^ Z^(x)£r  with  Y^Jx)  = T(x)Z^(x). 

r=  0 

In  particular,  therefore,  (10.4)  implies 


(10. 12) 


2q(x)  = e 


This  must  be  a solution  of  the  equation 


(10.13) 


x"k  df  = B0(x)z0 


Therefore, 


(10. 14) 


~ -L-  o -1  -k 

B0(x)  = x dT  Z0  = x Q,(x) 


In  explicit  form,  if  Q(x)  4 0,  then 


(10. IS) 


Q(x)  * 2i  QjX  ’ Qm  * ° ’ 
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so  that 


(10.16) 


- 171  . , . , . m 

R \ m -k-m-1  y ^ m~ 

- L jC^X  = ~X  /.  jQ  x 


j = l 


L 

j = l 


jQ; 


unless  Q(x)  = 0,  when  BQ(x)  = 0. 


Returning  to  (10.11),  one  sees  that  by  virtue  of  (10.8)  the  matrices 


(10.17) 


Br(x)x 


VV"r 


are  bounded  in  0 < |x|  <Xg,  I arg  x|  < 6^  although  not  necessarily 


holomorphic  at  x = 0.  Now,  set 
(10. 18) 

The  differential  equation  (10.9)  can  then  be  written 


B(x,e)  = B(x,  e)xk+m_1,  Br(x)  = Br(x)xk+m_1 


(10.19)  xm+1  ^ = B(x,  e)z  , 

and  the  B^(x)  satisfy  an  order  of  magnitude  inequality  of  the  form: 

pr  < 

(10.  20)  B^(x)x  is  bounded  for  0 < |x  | < xQ,  larg  x I < ; 

(10.21)  < r/«  , 0 < < °°  . 

* 

The  proof  that  B(x,  e)  i d , i.e.,  that 


(10. 22) 


B(x,  e)  ~ \,  B (x)er 


r=  0 


in  the  sense  of  Definition  2.1,  resembles  so  much  analogous  previous 
arguments,  e.g.  in  section  7,  that  a repetition  is  not  necessary. 

Thus  the  following  lemma  has  been  obtained. 
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Lemma  10 . 2.  The  transformation  (10-7)  (where  T(x)  is  defined  by  ( 10 ■ 6) 
and  (10 . 4)f  takes  the  differential  equation  (10.1)  (which  is  supposed  to 
have  Property  8)  into 


(10. 23) 


where  B(x,  e)  * & • 
( 10.  22)  for  B(x,  e) 

(10. 24) 


m+1  dz  , 

X dx  = B(X’ C'Z  ’ 


The  leading  term  B (x)  in  the  asymptotic  series 
is  either  zero,  and  then  m = 0.  or 


m 


Bq(x)  = - Yj  SQjx 
j = l 


m-j 


(with  Q.  as  in  (10. 15).  ) 


00 

\ 


The  coefficients  Z (x)  of  the  series  Z (x)e  which  solves 

r r = 0 r 

(10. 23)  formally  must  be  solutions  of  the  differential  equations 


(10. 25) 


, dZ 
m+1  r 


dx=B0(x)Zr+  £ B (x)Zv(x),  r = 0,1,... 
^+v=r  K 
^>0 


which  is  the  analog  of  (10.2). 

The  asymptotic  calculation  of  the  Z^x)  is  based  on  the  lemma  below. 

. . 2 , .m 

Lemma  10.  3.  Let  q(t)  be  a scalar  polynomial  q(t)  = q^  + q^t  +---+qmt  , 

m > 1.  a * 0 . and  let  q(t)  be  a scalar  function  holomorphic  and 
— ’ m 

bounded  in  the  region 

S = {t | tQ  < It  I,  larg(t  - 0Q) I < 6Q}  . 

Then  the  function 

f(t)  = eq^g(t)th  (h  an  integer) 


I 
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has  an  antiderivative  F(t)  of  the  form 

F(t)  = eq(t)G(t)th"m+1 

with  G(t)  holomorphic  and  bounded  in  some  subsector 

jjt  I 

S = {tjtj  < It  I , I arg( t - 0Q)  I < 6^ } 
of  S.  (G(t)  has  nothing  to  do  with  the  matrix  G in  (10.4).) 


Proof.  Choose  t^  > t^  so  large  that  the  mapping  s = q(t)  has  a 
holomorphic  inverse  in  the  subsector  of  S with  It  I > t^.  Define  F(t) 
by  the  integral 


F(t)  = J eq^g(T)-rhdT  . 
v(t) 

The  path  yU)>  which  must  end  at  t is  to  be  suitably  chosen  later. 

Since  t = q ^ms^m(l  + 0(s  *^m)),  as  s — »,  the  transformation 
m 

q(T)  = a changes  the  integral  into 

F(s)  = F(q(t))  = F(t)  = / eag(<T)«x(h"m+1)/mda  . 

\(s) 

A 

y(s)  is,  of  course,  the  image  of  the  path  y(t)  with  g(a)  bounded 
and  holomorphic  in  the  image  S = q(S)  of  S.  When  written  as 

(10.26)  F(s)e”ff  = / ea_Sg(1T)<r(h'm+1)/md<T  , 

v(s) 

the  integral  takes  on  a form  dealt  with  in  Lemma  14.  2 of  [18),  where  the 
setting  is  more  general,  except  that  the  exponent  that  corresponds  to 
the  exponent  of  a in  the  integral  above  is  assumed  to  be  non-positive. 
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This  assumption  is,  however,  nowhere  necessary  in  the  proof.  Instead 


of  repeating  the  rather  lengthy  arguments  of  that  proof  it  will  suffice 
to  state  the  conclusion  which  is  that  the  integral  in  (10  .26)  is  0(s^ 
as  s - or  m q(s  ),  if  6^  is  taken  small  enough.  This  completes 
the  proof  of  Lemma  10.  3. 

Lemma  10 . 4.  The  differential  equations  (B.29)  possess  solutions  of  the  form 

(10.27)  Z (x)  = Z.(x)eQ(x)  , 

/N.  v “Cl 

(10.28)  Z (x)  = Z(x)x  r, 

where  Z^(x)  is  holomorphic  and  bounded  in  0 < |x  I < xQ,  !arg  x I < 6 , and 

(10.29)  <r  < t/k 

for  some  number  k in  0 < « < <* . 

Proof.  Formula  (10.12)  establishes  the  lemma  for  r = 0.  For  r > 0, 
replace  (10.29)  by  the  equivalent  integral  representation  of  the  Z^ 

x . 

(10.  30)  Zf  = J ZQ(X)Z“  (i)  £ b (4)Z  l&)rm'  • 

p+v=r  ^ 

p>0 

The  integral  sign  is  to  be  interpreted  as  indicating  an  as  yet  unspecified 
antiderivative  taken  at  § = x.  Assume  the  lemma  to  have  been  proved 


for  Zq,  Zj,  . . . , Z ^ By  making  use  of  this  hypothesis,  as  well  as 
formulas  (10.20),  (10.21)  and  (10.12)  one  finds,  after  multiplication  with 
e on  both  sides  of  the  equation,  that  (10.  30)  can  be  written  in  the  form 

(10.  31)  Zf  = / eQ(x)-Q(^Hr(|)eQ(^"Q(x)i  rd£  . 
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A short  calculation  shows  that  if  h satisfies  the  inequality 

r 

(10.  32)  h < 1 

r k 

where  k has  been  taken  so  small  that 

k £ Kj/l  1 + (m  + 1)kj  ] , 

with  Kj  as  in  (10.  21),  then  H^(x)  is  bounded,  as  x — 0 in  larg  x I < 6^. 

The  right  member  of  (10.  31)  is  a matrix  of  n^  scalar  integrals 
of  the  form 


(10.  33) 


/ VX)"V£) 

J e 


h.k(e)rrAde,  j,k  = l,  ...,n. 


The  q (x)  are  the  differences  of  the  polynomials  in  1/x  that  appear 
J K 

as  the  diagonal  entries  of  Q(x).  Their  degrees  may  therefore  be  any 

integer  between  0 and  n.  If  the  degree  of  q.,(x)  is  positive, 

J K 

Lemma  10.  3 can  be  applied,  after  the  transformation  £ = 1/T-  The 

2-t/k 

result,  in  terms  of  x is  then,  that  the  integral  is  0(x  ),  as 

x -►  0.  If  q (x)  is  constant,  an  elementary  argument  shows  that  the 
I-t/k 

integral  is  0(x  ),  except  if  r/x  = 1,  which  can  be  avoided,  by 

~ -t/k 

not  choosing  « as  an  integer.  Thus,  the  relation  Z^(x)  = 0(x  ), 

as  x -*  0 in  larg  x I < is  always  true,  and  the  lemma  is  proved. 
The  matrices 


(10.  34) 


Yr(x)  = T(x)Zr(x),  r = 0,1,... 


solve  the  differential  equation  (10.2).  It  remains  to  be  shown  that  the 
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formal  series  solution 


or 

^ Y (x)er  represents  asymptotically  a solution 
r=  0 r 

of  (10.1).  It  is  convenient  to  prove  first  the  analogous  statement  for 

oo 

\ ' r 

the  senes  Z (x)e  with  respect  to  the  equation  (10.23).  The 

r=  0 r 

return  to  (10.1)  is  then  a simple  matter. 

With  the  notation  of  Lemma  10. 4 let  Z be  a function  in  Q 
with  the  asymptotic  expansion 


(10.  35)  Z(x,  e)~  \ Z (x)er,  as  e -*  0+  , 

r=0  r 


for  xi  £ . Such  a function  exists  according  to  Lemma  2.4,  which 


remains  literally  valid  in  O instead  of  Ct. 
Since 


00 

y 


r=  0 

(10.23),  it  follows  that 


Z^(x)e^^X^er  is,  by  construction,  a formal  solution  of 


(10.  36) 
where  *(xfe)~0,  i.e., 

(10.  37) 


(xm+1  - B(x,  e))(Z(x,  e)eQ^)  = *(x,  e)eQlx^  , 


*(x,  e)  < cN(  |x  I 1/ke)N+1  , 


for  x < & , 0 < e < and  for  all  N > 0.  Therefore,  if  (10.23)  is 
e 0 — 

converted  into  a differential  equation  for  the  matrix-valued  function  V 
by  inserting 

(10.  38)  z = (V  + Z(x,  e))e^X^ 

into  it,  one  gets,  with  the  help  of  (10. 15),  ( 10 . 24)  and  (10 . 36),  the 
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differential  equation 


(10.  39)  xm+1  ^ = BqV  - VBq  + (B  - BQ)V  - *(x,  e)  , 

and  the  task  is  now  to  show  that  it  has  a solution  that  is  asymptotic 
to  zero.  To  achieve  this,  one  observes  that,  according  to  a version 
of  the  variation  of  parameters  formula  derived  in  [ 18  ] , § 30.1,  every  solution 
of  the  integral  equation 

(10.40)  V(x,  e)  = / eQ{x]'QU)[(B(e,e)-B  (4)]V(4,e)-*(x,  E)]4“m_1eQ(U"y(x)de 

Hx) 

also  solves  (10.39).  Here,  formulas  (10.15)  and  (10.24)  have  again 
been  used. 

The  estimation  of  the  integral  operator  in  the  right  member  of  (10.  40) 
begins  with  the  observation  that 

Ib(|,  e)  - dQU)l  < c0UrlAe  , 

for  i t 0 < e < Eq,  in  consequence  of  (2.6).  Now,  replace  x 

* * 

by  a smaller  number  k , 0 < k < k,  such  that 

* < * 

K 1 + x(m  + 1)  ’ 

* 

and  write  f)  for  the  subset  of  obtained  by  this  change  in  the 

E E 

definition  of  fl  . Then 
£ 

(10.41)  Ib(4,e)  - B0(4)l  III""1'1  < c0tJ)/K+m+1E1“'<  (1/,<+m+1) 

* 

for  x < ft  , 0 < e < eq,  and  the  right  member  can  be  made  as  small  as 
one  pleases  by  reducing  eQ  or  tQ. 


-66- 


Let  the  q , (x),  j,  k = 1,  2,  . . . , n be  defined  as  in  (10.  33),  then 

the  set  of  n"  paths  of  integration,  each  ending  at  £ - x must  be  such 

that  in  each  scalar  integral  in  (10.40)  the  function  Req  (|)  is  non- 

1 k 

increasing.  The  proof  that  this  can  be  done  for  x t f.^,  provided  tQ,  6Q, 
e are  chosen  small  enough  is  a standard  argument  that  will  be 

omitted . 


A solution  V of  (10.  40)  which  is  asymptotic  to  zero  can  now  be 
constructed  by  the  usual  Picard  iteration:  Set  V = 0,  then 


V (x,e)  = - / 

1 r(x) 


The  path  of  integration  has  finite  length,  the  exponential  factors  are 
bounded  for  the  appropriate  choice  of  the  matrix  of  paths  P'W,  and 
(10 . 37)  implies 


I./*  x*-m-l|  * , | | -1/ k ^N+l 

l*(£,  e)£  I < CN(  lx  I e) 


* 

for  x i it  , 0 < e < e„  and  all  N > 0.  (C.r  is  a constant).  Hence,  also 
e o ~ n 


* 


(10. 42) 


v^x,  e)  I < CN(  lx  I 1/k  e)N+i 


with  some  constant  C^. 

The  details  of  the  proof  from  here  on  - whether  based  on  the 
Contraction  Theorem  in  Banach  Space  or  explicitly  on  the  Picard  iteration 
argument  - are  analogous  to  those  found,  with  variations,  in  almost 
every  paper  on  the  asymptotic  nature  of  formal  series  solutions  tor 
differential  equations.  It  suffices  to  point  out  that  the  possibility  of 
making  the  right  member  of  (10.41)  small  is  the  basis  of  the  existence  and 
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uniqueness  proof  for  the  solution  of  the  integral  equation  (10.40),  while 
(10.  42)  implies  that  this  solution  is  small  of  the  same  order,  i.  e. , that 

(1°-  43)  V(x,  e)~0  for  x € jo  0<e<e  , 

e - 0 ’ 

in  the  sense  of  Definition  (2.6). 

1 inally,  the  last  result  must  be  translated  into  the  equivalent 
statement  for  the  differential  equation  (10.1)  by  the  transformation  (10.7), 
replacement  of  x by  x in  the  notation,  and  return  to  x by  the 
inverse  of  the  change  of  variables  (10.  5).  The  proof  of  the  theorem  below 
is  now  complete. 

Theorem  10.1.  The  differential  equation  (10.1) : 

-k  dy  .. 

X dx  = A(X’  e)y’  A(X’  e)  4 

k an  integer  - positive,  negative  or  zero  - has  a fundamental 
matrix  solution  of  the  form 

(10.43)  Y(x,  e)  = Y(x,  e)xGeQ(x) 

with  the  following  properties. 

(i)  G and  Q(x)  are  the  matrices  defined  in  Lemma  10.1; 

(ii)  Y(x,  e)  = YQ(x)Y(x,e)  , where  YQ(x)  was  described  in 

*>*  * 

Lemma  10.1.  and  Y(x,e)  t C7  , 

lim  Y(x,  e)  = I . 
e -*0  + 

~ Q jJj  jjt 

Remark.  The  product  Y(x,  e)x  is,  of  course,  in  class  a . Only  this 
weaker  fact  is  needed  for  the  Main  Theorem  stated  in  the  Introduction. 
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§ I 1 . The  Main  Theorem. 


The  several  successive  transformations  of  the  original  differential 
equation  (1.1)  can  be  assembled  into  one  transformation 

(11.1)  y = T(x,  e)z,  x - x1",  e = e*  . 

T(x,  e)  is  a product  of  matrices  each  corresponding  to  a transformation 
of  Type  I,  II,  III  or  VI.  Whenever  one  of  these  transformations  is  applied 
to  one  of  the  blocks  of  the  block  diagonalization  reached  at  that  stage 
it  can  be  extended  to  a transformation  of  the  whole  space  which  leaves  the 
remaining  components  unchanged.  Observe  that  the  transformations  of 
Type  VI  are  scalar  in  each  of  the  blocks  of  all  subsequent  transformations. 
Thus,  these  transformations  can  all  be  moved  to  the  rightmost  positions. 
Therefore,  T(x,  e)  in  (11.1)  has  the  form 

(11.2)  T(x,  y)  = T(x,  e)exp{e  KQJ(x,  e)} 

with  the  same  structural  properties  as  formula  (1.  3). 

The  differential  equation  for  z resulting  from  transformation  (11.1) 
consists  of  scalar  uncoupled  equations  with  solutions  of  the  form 
described  in  Theorem  9-1  and,  jiossibly,  of  systems  without  a power  of  e 
multiplying  the  derivative.  These  systems  were  solved  in  section  10.  The 
solutions  so  obtained  can  be  assembled  into  a matrix  Z(x,  e)  of  the  same 
form  as  (11.2).  Moreover,  the  exponential  factor  in  (11.2)  commutes  with 
Z,  since  the  diagonal  matrix  Q (x,  e)  is  scalar  in  each  diagonal  block 
of  Z.  Therefore  the  exponential  factors  in  the  product  Y = TZ  can  be 
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combined,  which  completes  the  proof  of  the  Main  Theorem  as  stated  in 
the  Introduction. 

The  equations  of  the  form  (1.1)  for  which  the  outer  solution,  as  described 
in  Theorem  1.  1 has  been  explicitly  connected  with  an  interior  or  intermediate 
solution  (see,  e.g.  [12],  [13],  [14],  [ 19  | , [22])  all  possess,  to  my 
knowledge,  the  special  property  that 


(11.  3)  le  HQ(x,  e ) I < K 

for  some  constant  K and  for  certain  values  of  x and  e that  lie  in 

the  domain  . 0 < e < in  which  the  asymptotic  results  of  Theorem  1.1 
e — 0 

* 

are  valid.  This  is  by  no  means  always  the  case,  even  if  Q « (7  , as 
will  be  illustrated  by  the  example  below.  That  example  also  shows  that 
the  methods  of  this  paper,  while  constructional  in  principle,  involve 
long  computations,  even  for  simple  equations. 

Examples . 

The  equation 


(11.4) 


2 dy 
dx 


0 1 0 
0 0 1 


e 0 x 

is  of  the  type  (1.1)  with  h = 2 


m 


and 


m a positive  integer  , 


f°  1 0 ) 

1 0 0 0 \ 

(11.  S)  A(x,e)  = A0(x)  + Aj(x)e  = 

0 0 1 
(o  0 xmi 

+ 

— o 
o o 
o o 

while  the  coefficient  matrix  is  a simple  companion  matrix,  the  turning 
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point  at  x - 0 is  relatively  complicated.  Two  eigenvalues  of  A^fx) 
are  identically  zero,  and  the  third  eigenvalue  becomes  equal  to  the 
other  two  at  x - 0.  The  elementary  divisors  of  A^( x)  have  degrees 
2 and  1 for  x 0,  but  for  x = 0 there  is  only  one  elementary  divisor, 
of  degree  3. 

As  ( 0)  is  already  in  Jordan  form  and  nilpotent,  the  decomposition 
of  the  system  begins  with  the  shearing  described  in  § 7 . The  figure 
below  shows  the  Puiseux  polygon  for  equation  (11.4). 


v 


The  appropriate  shearing  transformation  is 

(11.6) 


,,  m 2m,  ~ 
y = diag(l,x  ,x  )y  , 


and  the  transformed  differential  equation  becomes 


/ii  2 -m  dy 
(11.  7)  e x , = 
dx 
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Next,  the  transformation 
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^ -w. 
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y - T y with  T = 
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(11.8) 


block-diagonalizes  the  new  lead  matrix.  One  gets 


(11.9) 

dx 
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The  formal  block  diagonalization  of  the  whole  coefficient  matrix  by  the 
technique  of  § 4 introduces  a non-terminating  series  in  powers  of  e.  The 
calculations,  while  elementary,  are  tedious.  One  finds  that  a transforma- 
tion of  the  form 


(11.10) 

with 

(11.11) 


y = P(x, e)z 


P(x, e)  = I + 


1 0 0-1^ 


0 0-2 

-1  -1  0 


-3m  . -6m  2. 

x e + 0(x  e ) 


takes  the  differential  equation  (11.9)  into  the  block-diagonal  form 


/it  i , t ^ — m d z 

(11. 12)  e x — = 
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which  splits  into  the  two  uncoupled  problems 


2 -m  du 
: x T~ 
dx 


( 

jo  l) 

1-1  o\ 
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1 + 

l-i  o j 

r e+i 

\l  1 - mx 


5m-  1 


-6m  2 . -9m  3 

x e +0(x  e 


......  2 -m  dv  , -3m  , , . _ 5m-l,  -6m  2 , -9m  3n 

(11.  14)  e x — = { 1 + x e + (-1  - 2mx  )x  e + 0(x  e )}v  . 

T 

Here  z = (u,  v)  . 

Tirst,  the  analysis  of  (11.13)  will  be  continued.  As  the  reduction 
to  Arnold's  form  according  to  § 6 is  not  essential  here,  it  will  be  skipped. 
Another  shearing  is  then  in  order.  Below,  the  points  of  the  Puiseux 
diagram  which  can  be  recognized  from  the  partial  sum  in  (11.13)  are 
plotted.  From  the  torm  of  the  series  it  is  clear  that  the  remaining  points 
are  irrelevant  for  the  exponents  of  this  shearing. 


v 
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The  indicated  shearing  transformation  for  (11.13)  is 


(11.  IS) 


,,  1/2  -3m/2.~ 

u = diag(l,  £x  ) u . 


It  produces  the  new  differential  equation 

„ ...  2 2d=  [1°  ‘\l‘  °\  'I1"  1/2  M - 3m  1°  M 3/2 

(11. 16)  e x — =(  + x e + x e+  x e 

dX  U-l  0/  \0  1/  \0  1 / \0  -m / 

-0m  2.  1 ~ 

+ 0(x  e ) ) u . 

The  re-adjustment  to  integral  exponents  will  be  bypassed  as  inessential, 
since  (11. 16)  can  be  diagonalized  directly  by  the  method  of  § 4.  The 
result  of  this  calculation  is  - details  are  omitted  - that  a transformation 
of  the  form 

(11. 17)  u = | J + 0(x'3me)  u 

changes  (11.17)  into  a diagonalized  system  of  the  form 


(11.18)  e 3/2xm/2  J1  = 

dx 
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Now  the  diagonal  system  (11.18)  and  the  scalar  equation  (11.14) 
must  be  integrated  by  the  elementary  procedure  of  §9  and,  finally,  the 
several  transformations  performed  in  this  section  can  be  combined  to 
yield  the  final  result.  Neither  the  details  of  the  final  formula  nor  the 
calculations  have  much  intrinsic  interest,  except  for  the  structure  of 
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e — 0,  in  some  parts  of 

The  properties  of  the  differential  equation  (11.4)  in  the  shrinking 
disk  around  x = 0 are  very  complex.  Details  of  a general  analysis  - 
though  not  a solution  of  the  problems  that  arise  - are  given  in  [10] 
and  [6],  ( 7 ] . 

The  natural  tool  is  a "stretching"  transformation 
(11. 19)  x = 4ep,  p > 0 . 

It  changes  (11.4)  into 
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0 ] 0 
0 0 I 

e 0 4 

If  p < l/3m  any  finite  disk  in  the  £-plane  overlaps  with  at  least 

for  small  enough  e.  However,  as  one  can  check  directly,  the  theory 

of  this  paper,  when  applied  to  (11.  20)  does  not  extend  the  validity  of 

the  asymptotic  analysis  to  points  closer  to  the  origin  than  0(xe  ' ). 

For  p - l/3m,  the  shearing  indicated  by  the  theory  is 
, 1/3  2/3, 

y = diag(l,  e , e )z  and  leads  to 

0 1 0 
0 0 1 
1 0 

Now,  the  origin  is  no  longer  a turning  point,  but  there  are  3m  isolated 
points  in  the  4-plane  where  some  eigenvalues  of  the  coefficient  matrix 
coalesce.  Thus,  the  local  problem  in  the  x-plane  has  become  a very 
complex  global  problem  in  the  £-plane. 

The  occurrence  of  such  secondary  turning  points  after  the  stretching 
is  typical,  whenever  the  original  Puiseux  polygon  has  more  than  one  side. 
The  paper  [11]  by  M.  Nakano  and  T.  Nishimoto  illustrates  well  these 
complications. 


( 5m-2)/3  dz 
E d£  Z 
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